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SLTMMARY 


The general design method developed in this report is for steady, three-dimension- 
al, potential, incompressible or subsonic -compressible flow. In this design method, 
the flow field, including the shape of its boundary, is determined for arbitrarily speci- 
fied, continuous distributions of velocity as a function of arc length along the boundary 
streamlines. The method applies to the design of both internal and external flow fields, 
including, in both cases, fields with planar symmetry. These designs result from the 
finite-difference solution of a governing, partial-differential equation for the distribu- 
tion of velocity in transformed space, the coordinates of which are the velocity poten- 
tial and two stream functions. The analytic problems associated with stagnation points, 
closure of bodies in external flow fields, and prediction of turning angles in three- 
dimensional ducts are discussed, but not treated in detail. 

This three-dimensional design method applies to simple ducts, branched ducts, 
annular ducts, and ducts with centerbodies (including the nose and/or tail region) and to 
ducted or nonducted bodies in infinite, uniform flow fields. 


INTRODUCTION 

In the theory of multidimensional, potential flow there are two kinds of boundary- 
value problems: (1) the direct problem, in which the distribution of velocity is deter- 
mined for a prescribed shape of the boundary and (2) the inverse problem, in which the 
shape of the boundary is determined for a prescribed distribution of velocity along it. 
The direct problem is an analysis problem; the inverse problem is a design problem. 
This report is concerned with the inverse, or design, problem for steady, three-dimen- 
sional, subsonic, potential flow through ducts, or around bodies, with prescribed veloc- 
ities as a function of arc length in the direction of flow along the boundary of the field. 

The design of flow fields with satisfactory velocities along the boundary is impor- 
tant for the following reasons: 

(1) Boundary-layer separation losses can be avoided by prescribed velocity distri- 
butions in the direction of flow, along the material surfaces of the boundary, that do not 
decrease too rapidly. 

(2) Shock losses in compressible flow, and cavitation in incompressible flow, can 
be avoided by prescribed velocities that do not exceed certain maximum values dictated 
by these phenomena. 

(3) For compressible flow in ducts, the desired flow rate can be assured by pre- 
scribed velocities that do not result in premature choked flow. 


However, the first objective of fluid dynamic design is to determine the shape of 
the boundary of the flow field for which losses are minimum. For both incompressible 
and shock-free compressible flow, fluid losses originate at the material surfaces along 
the boundary of the flow field, and the magnitude of these losses depends on the velocity 
distribution along these surfaces. The characteristics of a desirable velocity distribu- 
tion are relatively well-known from boundary-layer theory. 

Solutions of either the design or analysis type of boundary-value problem can be 
used to obtain a satisfactory velocity distribution. In the inverse, design method, an 
acceptable velocity distribution is specified, and the resulting shape of the flow field 
determined. If for some reason the shape is unsatisfactory, the prescribed velocity 
distribution can be modified, and the solution repeated until an acceptable shape is 
achieved. In the direct, analysis method, on the other hand, an acceptable shape is 
specified, and the resulting velocity distribution determined. If this velocity distribu- 
tion is unsatisfactory, the shape can be modified until an acceptable distribution is 
achieved. 

With either method, the final design is usually arrived at only after iterating on the 
specified input. The inverse, design method, however, has some advantage, because 
it starts with an acceptable velocity distribution and because, in some cases, it arrives 
at boundary shapes that could not otherwise be envisioned by the designer. 

The previous discussion applies both to internal flow fields, in which the boundary 
surrounds the flow, as in ducts, and to external flow fields, in which the flow surrounds 
the bounding surface, as on bodies . 

The general design method developed in this report is for three-dimensional, poten- 
tial, incompressible or subsonic -compressible flow. The method applies to the design 
of both internal and external flow fields, including, in both cases, fields with planar 
symmetry. For plane, two-dimensional flow, design methods have been developed in the 
past for internal flow fields (e.g., ref. 1), and numerous methods have also been de- 
veloped for external flow fields. For axisymmetric, two-dimensional flow, a design 
method has been developed recently for internal flow fields (ref. 2), but no similar de- 
sign method appears to have been attempted for external flow. Until now, no general 
three-dimensional design method has existed for either internal or external flow fields. 

A major difficulty faced by all the design methods just mentioned is that they are 
boundary-value problems in which the velocity distribution is specified along physical 
boundaries, the shapes of which are not known until the problem is solved. In the two- 
dimensional cases, this difficulty has been avoided by solving the boundary-value prob- 
lem in the transformed plane of the velocity potential q? and the stream function In 
the present three-dimensional case, the difficulty is avoided by solving the problem in 
(jci, 77-space, where 77 is a second stream function associated with continuity in three- 
dimensional flow (ref. 3). 
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Because of its general nature, the method of this report is not likely to be used for 
the design of flow fields with special physical characteristics, such as the hub and 
shroud surfaces of revolution in turbomachines. These particular surfaces, for exam- 
ple, should be more easily designed by using the method developed in reference 2 for 
axisymmetric flow. However, as is demonstrated later in this report, the design 
methods for plane, two-dimensional flow (e.g., ref. 1) and for axisymmetric, two- 
dimensional flow (ref. 2) are special cases of the general, governing differential equa- 
tion developed in this report. 

Finally, the analyses of this first report do not treat problems associated with the 
following three important boundary conditions: (1) stagnation points, (2) body closure, 
and (3) duct turning. Thus, in this analysis, any region in the vicinity of what would 
normally be a stagnation point is cusped; body closure can be assured only when a plane 
of symmetry exists that is normal to the direction of undisturbed flow (otherwise the 
body will likely have a tail of essentially constant cross section that trails downstream 
to infinity); and the desired duct turning can be achieved only by trial -and -error itera- 
tion of the specified velocity distribution on the boundary of the flow field. 


SYMBOLS 


A 

a 

B 

C 


All quantities are nondimens ional unless otherwise specified. 

local continuity parameter, eqs . (10a) and 0.0g) 

upstream boundary area, 1.0 in nondimensional form 

distance between adjacent nodal points of finite-difference star, 
fig. 10 

local continuity parameter, eqs. (10b) and (lOe) 
constant 



c 


3’ ^4 


E,F,G 


e 


coefficients in governing finite-difference eq. (28e), eqs. (28f) 

coefficients given by eqs. (21f), (22f), and (23f) 
local speed of sound, m/sec, eq. (3b) 

coefficients of first fundamental quadratic form of surface geometry, 
eqs. (A3) 

unit tangent vector 


unit vector in direction of q along streamlines, which are intersec- 
tions of ip and rj stream surfaces, fig. 3 and eq. (7a) 
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H 

h 

i,j,k 


k 

M 

m 


n 


n 


P P 
1 ’ ^ 2 ’ 

. . .,P 


12 


PU 

Q 

q 

q 

R 

R 

(R 


r 

r 

T G * 

s 

T 


u,v 


unit vector tangent to intersection of 17 stream surface and tp poten- 
tial surface, fig. 3 and eq. (7c) 

unit vector tangent to intersection of ip stream surface and (p poten- 
tial surface, fig. 3 and eq. f7e) 

(8r/9u) X (aF/8v) • n, eq. (B2) 

enthalpy, JAg 

unit vectors in x-, y-, and z-directions, respectively 

total curvatures of 77 and tp stream surfaces, respectively, eqs. 

(17 f) and (17 g) (also (C5) and (C6)) 

nondime ns ional parameter, eq. (4e) 

local Mach number, eq. (3a) 

path length in direction of e^ along intersection of ip stream surface 
and (p potential surface, fig. 3 

path length in direction of Fg along intersection of 77 stream surface 
and cp potential surface, fig. 3 

unit vector normal to surface 

parameters given by eqs. (2lg), (22g), and (23g) 

path length along perimeter of upstream boundary, fig. 11 
In q 

velocity 

velocity vector, e^q, fig. 3 and eq. (lb) 
gas constant, J/^g)(K) 
defined by eqs. (C2) and (C4) 
residual error, eq. (28e) 
polar coordinate 

position vector in physical space 

polar coordinates on upstream boundary, fig, 12 and eq. (30g) 
path length along streamline in direction of or entropy, J/(kg)(K) 
absolute temperature, K 
curvilinear coordinates 
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u, v,w 


V 

x,y,z 


a,/3,7 


y 

€ 

V 

e 


0 ’ 

p 

<p 

(Pi^yT] 


V 

X 


Subscripts; 

C 

D 


max 

min 

o 

S 

u 

X 


components of velocity q in directions of i, j, and k, respectively 
volume 

Cartesian coordinates in physical space 
dummy variable 

angles of direction cosines in x, y, z-space, fig. 1 
ratio of specific heats; or geodesic curvature 
error 

stream function, eqs. (6c) 

angle with which ip and 17 stream surfaces intersect on potential sur- 
face, fig. 3 and eqs. (8b) and (8c) 

polar coordinate 

local density of fluid, eq. (4d) 

velocity potential, eqs. (5d) 

curvilinear coordinates in physical x,y, z-space or coordinates of 
transformed 77-space 

stream function, eqs. (6c) 

vector differential operator, eqs. (5b) and (12a) 
vector differential operator on surface 
cross-product operator of two vectors 
dot-product operator of two vectors 

cut line or surface cut 
downstream boundary 
maximum 
minimum 

total, or stagnation, condition 
surface 

upstream boundary 

any grid point on a potential surface in (^, ^, 77-space 
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0 , 1 ,. . .,6 grid points in finite-difference star, fig. 10 

1, 2, 3 variables or components of variables associated with directions of 

e^, eg, and respectively 

+, - positive and negative directions, respectively, from central grid point 

in relaxation star, eq. (28a) 

* dimensional quantity 


THEORY OF DESIGN METHOD 

The theory of the design method is developed in this section. After preliminary 
considerations, the governing, second-order partial-differential equation for distribu- 
tion of velocity in transformed ^,^,i 7 -space is developed; additional equations, making 
use of results from the solution of the governing equation, are then developed for con- 
struction of the corresponding flow field in physical x,y, z-space. 


Preliminary Considerations 

The purpose of this design method is to determine the boundary shape of the poten- 
tial flow field, in three-dimensional x,y, z-space, for which field the distribution of 
velocity everywhere on the boundary is as specified. This distribution of velocity is 
specified by q(s), where s is distance along streamlines on the boundary . But, for 
potential flow, the velocity potential cp (s) varies along each streamline according to the 
well-known relation (see eq. (5f)) 

q(s)ds 

As a result, the prescribed distributions of velocity q(s) on the physical boundary, of 
unknown shape, become prescribed distributions of velocity q(cp) along the known cur- 
vilinear coordinate <p . If two additional curvilinear coordinates ^^,y, z)and y, z), 
related to the streamlines, can be determined, then the shape of the flow field in trans- 
formed cp y jp, T)- space is known, as is the prescribed velocity q(s,j/),? 7 ) = q(<^, on its 
boundary. The distribution of q throughout the flow field in cp, j/), 77 - space can then be 
determined by a governing, partial-differential equation, and, from this distribution, 
the coordinates of the flow field, including the physical shape of boundary, can be com- 
puted in x,y, z-space. In this subsection, the general curvilinear ^,t]-coordinates in 
x,y, z-space, which become Cartesian coordinates in 77 -space, are identified. 
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Problem statement . - Given the shape of the upstream boundary surface and the 
distribution of the velocity q as a function of the distance s along streamlines on the 
surface of a three-dimensional, potential flow field, find the shape of the field. 

Assumptions . - The flow is assumed to be steady and irrotational . The fluid is in- 
compressible or subsonic -compressible. (Although, in regions of limited extent, 
supersonic velocities can probably be accommodated.) 

Flow field in physical x,y,z-space . - A steady-flow field in x,y, z-space is defined 
by the spatial distribution of the velocity vector q and of the associated physical prop- 
erties of the fluid. The local velocity vector q is given by its magnitude q and by a 
unit vector e^ in the direction of q. This direction is tangent to the streamline, and, 
therefore, parallel to the elemental vector distance ^j^ds along the streamline. As 
shown in figure 1, at any point x,y, z, the direction of and thus of e^ds and ^^^q, 

is given by the direction cosines: cos ce^, cos and cosy^. Thus, if i, j, and k 
are unit vectors in the direction of x, y, and z, respectively, 

e^ = i cos + j cos + k cos y 1 (la) 


so that 


q = e^q = iu + jv kw 


m 


where u, v, and w are the velocity components in the x-, y-, and z -directions, re- 
spectively. Thus 


u = q cos ^ 


V = q cos /3^ 
w = q cos y^ 




(ic) 


Likewise, for ds. 


from which it follows that 


■>1 

dx = ds cos 

dy = ds cos (3^ > 

dz = ds cos y ^ 

J 


(Id) 


2 

cos + cos 


+ cos\^ = 1 


(le) 
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Density p . - For an incompressible fluid, the mass density p is constant, it does 
not affect the configuration of the flow field, and it, therefore, does not enter into the 
theory of the design method. For pure substances in the vapor phase only, that is, a 
real gas or a mixture of real gases, the static mass density p is a function of any two 
independent thermodynamic properties, of which at least one must be affected by the 
velocity q. It is convenient to select entropy s and static enthalpy h for these two 
properties . The static enthalpy is then related to the constant stagnation enthalpy h^ 
by 


h = h 

® 2 


(2a) 


and, for the irrotational motion of this design theory, the entropy s is constant, (More 
precisely, in the absence of external forces, and neglecting viscosity and heat conduc- 
tion or radiation, Crocco’s equation (ref. 4) shows that the flow field is irrotational, 
and, therefore, potential, if the entropy s and the total enthalpy h^ are constant.) 

For perfect gases with constant entropy, the static mass density is related to the 
velocity and to the stagnation mass density p^ by 


P = Po'l 


-T, r.i 




yRT 


(2b) 


where the constants y, R, and are the ratio of specific heats, the gas constant, 
and the stagnation temperature, respectively. 

Mach number M . - The local Mach number M is the velocity q divided by the 
local speed of sound c. Thus 



where 


c2=(M 

\dp/^gentropic 

Thus, from equations (3a) and (3b), for real gases. 


(3a) 


(3b) 
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For a perfect gas, 

=yRT 

\dp /jg entropic 

so that, from equations (3a), (3b), and (3d), 


(3c) 


(3d) 



where the local static temperature T is related to the constant stagnation temperature 
To by 



The local Mach number M is of critical interest in application of the design 
method to compressible flow fields, because the maximum velocity in the field, which 
velocity always occurs on the boundary and is, therefore, easily controlled, should 
normally be subsonic. Small local regions of supersonic flow can usually be accommo- 
dated, but strictly speaking, the fundamental character of the flow is different; the gov- 
erning partial- differential equation changes from an elliptic to a hyperbolic type. 

Nondimens ional forms . - From this point on, it is convenient to treat all variables 
as nondimens ional. Thus linear dimensions, such as x, y, z, and s, are normalized 
by the square root of the upstream boundary area in x, y, z -space: 



in which the subscript ♦ indicates a dimensional quantity. The velocity q and its 
components u, v, and w are normalized by the uniform velocity qjr at the upstream 
boundary: 

q = , etc. (4b) 

‘’u,* 


9 


And the density p is normalized by the constant static density ^ at the upstream 
boundary: 


P = 


^U, 


(4c) 


Thus, for incompressible flow, the density is 1.0, and, for compressible flow of a per- 
fect gas, from equation (2b), the density becomes 


P- 

\ 1 - k 

where q is the nondimensional velocity, and 




(4d) 



(4e) 


Also, for a perfect gas, from equations (3e) and (3f), the Mach number M, expressed 
in terms of q, becomes 


M = 



1 - kq^ 


(4f) 


from which, at the upstream boundary where M is and q is 1.0, 


k = 


y - 1 
2 

2 




(4g) 


so that k is zero when is zero, that is, when the flow is incompressible. Equa- 

tion (4g) determines k for a specified value of M^. 

Velocity potential function cp y, z) . - For irrotational flow, 

V X q = 0 (5a) 

where the vector differential operator V, expressed in terms of x,y, z-coordinates, is 
given by 
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(5b) 


so that, from equation (lb). 


— "T 9 "7 8 7” 8 

V = i — + j — + k — 


8x 




8z 


8w _ ^ _ Q ^ 

8y 8z 


8u _ 8w _ Q 
8z 8x 




8v 8u 
8x 8y 


(5c) 


From equations (5c), it follows that a velocity potential function (^(x, y, z) exists in the 
flow field such that 


y (5d) 


j 

The existence of cp^,y, z) is confirmed by substitution of equations (5d) into equa- 
tions (5c). Surfaces of constant cp in x,y, z-space are potential surfaces. 

Because ^,y, z) is a continuously differentiable point function of x, y, and z 
within the boundary of the flow field, 

d(^ = dx + dy + dz (5e) 

8x 8y 8z 

or, from equations (Ic), (Id), and (le), for an elemental distance ds along a stream- 
line, 

■^ = q (5f) 

ds 


U = §£ 

8x 


V = §£ 


W 


— 

8z 


Equation (5f) is the result of irrotational flow. 
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Finally, from equations (5b) and (5d), 


v<p=T^+7^ + k^ 

9x 8y 9z 


= iu + jv + kw 


or, from equation (lb), 


Vcp =q = e^q (5g) 

Because Vcp, the vector gradient of (ps is everywhere normal to surfaces of constant 
(p, the velocity vector q, which from equation (5g) is equal to V(p, must also be nor- 
mal to the surfaces of constant </?, as shown in figure 2. 

Stream functions ip y, z) and 7? (x, y, z) . - For steady, compressible flow, the 
continuity condition is 


V • (pq) = 0 

so that, from equations (lb) and (5b), 

afcu) , 8(pv) , 3(pw) _ Q 
8x 8y 8z 

where, for incompressible flow, the density p is 1.0. From equation (6b), it follows 
that two stream functions, ^(s,y, z) and ?7 (x,y, z), exist in the flow field such that 
(ref. 3) 


(6a) 


(6b) 


dij) dn __ ^ ^ 
dy dz dz dy 


pv 


M M 

dz dx 9x 9z 




= M M- MM. 

9x ^ 9y 9x J 


(6c) 


The existence of ^(x,y, z) and rj pz,y,z) is confirmed by substitution of equations (6c) 
into equation (6b). Surfaces of constant ip and surfaces of constant 17 are stream sur- 
faces in X, y, z- space . 

Finally, from the cross product of the gradients of ip and 77, 


12 



I 


V!j>xVr]=(i^ + }^ + k x(i^ + j^+k 

\ax dy 9z/\ax dy dz/ 




_ 


i 

j 

k 




9x 

dy 

9z 




dx 

dy 

dz 




dy 9z 9z 9y 


)'■( 


dip do _ 
dz dx dx dz 


H 


^ M 

dx dy dy dx 


) 


or, from equations (lb), (5g), and (6c), 


Vzp X Vtj = i(pu) + j (pv) + k(pw) 


= pq 

= pe^q 

= pVcp 




(6d) 


From equation (6d), the gradients Vip and Vt], which are normal to the stream sur- 
faces of constant ^ and of constant t], respectively, are normal to q and Vcp. It 
follows directly, as shown in figure 3, that the gradients of ip and of rj lie on surfaces 
of constant cp and, therefore, that stream surfaces of constant ip or of constant 77 are 
normal to potential surfaces of constant cp . (Also see appendix A.) Because the 
stream surfaces of constant ip or of constant 77 are both normal to potential surfaces 
of constant cp, their intersections are also normal to surfaces of constant (p and are, 
therefore, parallel to the velocity vector q. Thus the intersections of stream surfaces 
of constant ip with stream surfaces of constant 77 are streamlines, and the values of 
ip and 77 along a streamline are constant. In a similar fashion, the intersections of 
adjacent stream surfaces of constant ip and ip + 6 ip with adjacent stream surfaces of 
constant 77 and 77 + 677 define stream tubes, as shown in figure 4. Thus fluid does not 
flow across stream surfaces of constant ip and 77; these stream surfaces are, in ef- 
fect, generated by fluid lines of connected particles on the upstream boundary that pass 
downstream through the flow field and thereby divide it. Because the selected (input) 
shape and orientation of these fluid lines at the upstream boundary are somewhat arbi- 
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trary, the shape of the resulting stream surfaces and the angle 6 (fig. 3 ) with which the 
stream surfaces of constant ip intersect stream surfaces of constant 17 (measured on 
potential surfaces of constant <p) vary with fluid line selection for the same boundary 
shape of flow field in x,y, z-space. 

Curvilinear cp , ^-coordinates in x,y, z-space . - The two families of stream sur- 
faces of constant ip and of constant rj together with the family of potential surfaces of 
constant (p form a curvilinear ^,97 -coordinate system for the flow field in x,y,z- 

space. On stream surfaces of constant ip (fig. 3 ), there is a coordinate system of 97 
and <p lines. For these 77, (^-surfaces of constant the 17- and (^-coordinate lines 
are orthogonal, as shown previously by equation (6d). (Also see appendix A.) Like- 
wise, on stream surfaces of constant 17 (fig. 3 ), there is a coordinate system of cp 
and jp lines. For these <^, ^-surfaces of constant 97, the (p- and j^-coordinate lines are 
again orthogonal, as shown previously by equation (6d). (Also see appendix A.) On the 
potential surfaces of constant <p (fig. 3 ), there is a coordinate system of \p and 77 
lines. For these 97-surfaces of constant however, the and 97 -coordinates are 
not generally orthogonal, but intersect at an angle 0(<^,^,97) that varies with location in 
the flow field. (See fig. 3 and appendix A.) Thus, in this one respect, the curvilinear 

97-coordinates are not orthogonal in x,y, z-space. 

Unit vectors e^, and e^. - At every point of the flow field in x,y, z-space, the 

unit vectors and e^ (fig. 3 ) are tangent to intersections of the stream surfaces 

and the potential surface. The unit vector e^^ is tangent to the intersections of ip and 
97 stream surfaces and points in the direction of increasing distance s and of increas- 
ing velocity potential q) along the intersections (which are streamlines). This unit 
vector is defined in terms of its direction cosines by (also see eq. (la)) 


e^^ = i cos + j cos / 3 ^ + k cos 
where, from equations (Id), 



C7a) 


( 7 b) 


in which ds (fig. 3 ) is elemental distance along the intersections, and dx, dy, and dz 
are corresponding elemental distances in the x, y, z-coordinate directions. 
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The unit vector e„ is tangent to the intersections of 17 stream surfaces and (p 
potential surfaces and points in the direction of increasing distance n and of increasing 
stream function jp along these intersections . This unit vector is defined in terms of 
its direction cosines by 


e 2 = i cos Q ?2 3 ^2 ^ ^ "^2 


(7c) 


where 


cos = 

dx 

2 

dn 

cos = 


2 

dn 

cos y = 

dz 

2 

dn 


C7d) 


in which dn (fig. 3) is elemental distance along the intersections. 

The unit vector is tangent to the intersections of (p potential surfaces and ip 
o 

stream surfaces and points in the direction of increasing distance m and of increasing 
stream function rj along these intersections . This unit vector is defined in terms of 
its direction cosines by 


where 


e„ = i cos tto + j cos j3o + k cos y 


cos OJg = 


dx 


cos (3, 


cos y q = 


dm 

= ^ 
dm 

dz 


dm 




C7e) 


(7f) 


in which dm (fig. 3) is elemental distance along the intersections. 

The unit vectors and are tangent to the curvilinear coordinate lines 

and pointed in the directions of increasing cp, ip, and 77 , respectively. 
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Angle 0 . - Because the 
surface, the unit vectors 


jp and 7] stream surfaces are normal to the cp potential 
and e^ are normal to the unit vector e^ . Thus 



(8 a) 


And, because the two stream surfaces intersect at an angle 6(cp,iptT)) (fig. 3) on the 
^,Tj "Potential surface, 


63 • 63 = cos 0((^,j/j,T7) (8b) 

From equations (7c), (7e), and (8b), 

cos 6{(p,i{),r]) = cos cos O3 + cos cos (3^ + cos cos 7 3 (8c) 

(Also see appendix A.) 

The angle 0 is a major parameter in the governing, partial-differential equation 
of the design method. From figure 4 (and fig. 6), it can be seen that the range of d 
values is limited to 

0 < 0 < 180® (8d) 

Otherwise the area of the stream tube would become negative. 

Gradients Vq?, V^, and Vtj . - The gradients of cp., ipy and r\ in x, y, z-space are 
differential invariants that can be expressed in terms of any convenient coordinate sys- 
tem in X, y, z-space, for example, in terms of the natural coordinates s, n, and m 
(fig. 3). The gradient V(p, which is normal to potential surfaces of constant cp, is in 
the direction of e^ and has the magnitude d(p/ds. Thus 

V(^ = e^ ^ (9a) 

^ ds 

which also follows from equations (5f) and (5g) . 

The gradients and V?7, which lie on potential surfaces (fig. 3) and are normal 
to stream surfaces of constant ^ and to stream surfaces of constant r\, respectively, 
have the directions e 3 x e^^^ and 6^X62, respectively, as shown in figure 5 . Also, 
from figure 5, the magnitudes of these gradients are (1/sin d)d^/6n and (1/sin 0)dr]/dm, 
respectively, where sin 0 appears because d^/dn and d77/dm are in the directions of 
Cg and 63, respectively. Thus 
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I 


Vj/) = e„ X ^ ^ 
sin 0 dn 


(9b) 


V77 = e X e — ^ 

sin 0 dm 


(9c) 


These gradients are required to transform the fundamental equations of flow from 
x,y, z~space to (^, ^,rj-space. 

From equations (9b) and (9c), 


Vip xVtj ^ X e^^ X x 62^ 


sin^0 


M JiL 

dn dm 


where, from figure 5, 


(eg X e^) X x eg) = 


sin 0 


so that, from equation (6d), 


e^pq 


- 1 M M. 

^ sin 0 dn dm 


or 


dip dr] = pq(dn dm sin 0) (9d) 

Equation (9d) is the result of continuity, where pq(dn dm sin 0) is the differential mass 
flow rate across the elemental area (dn dm sin 0) bounded by adjacent ip and t] lines of 
differential magnitude dip and drj, respectively, on the potential surface in figure 6. 
From equations (5f) and (9d), it follows that 

dy d’f’ d’l = pq2 (9e) 

dV 

where, from figure 7, the elemental volume dV, bounded by adjacent potential and 
stream surfaces of differential magnitude d<p, dip, and drj, is given by 

dV = ds (dn dm sin 0) (9f) 

Thus, from equation (9e), for given values of d^, d^, and dr], the smaller the elemen- 

2 

tal volume dV, the greater pq . 
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From figure 5, the cross products involving e^ in equations (9b) and (9c) can be 
expressed in terms of the unit vectors and e^ lying on potential surfaces. Thus 


es X = 


Cl X 62 = 


— — - e« cos 0^ 

sin 0 \ ^ / 

— i e n - e^ cos 0^ 
sin 0 \ / 


from which, equations (9b) and (9c) become 


V?7 


sin^0 




e - e„ cos 0\ ^ 
^ ^ / dn 


e^ - cos 0 


/ dm 


(9g) 

(9h) 


(9i) 

(9j) 


Parameters A and B . - In the development of the design method, it is convenient 
to introduce two parameters, A and B, which are defined by 


A = p — sin 0 (10a) 

d^ 


B = p — - sin 0 (10b) 

di7 

In equation (10a), and as shown in figure 8, dn sin 0 is the local thickness of the dif- 
ferential stream lamina bounded by the adjacent stream surfaces of constant ip and of 
constant ip + dip. The parameter A, which is this thickness multiplied by the local den- 
sity, is, therefore, a measure of the local mass-flow capacity of the stream lamina for 
any given value of the local velocity q. Likewise, from equation (10b) and figure 8, the 
parameter B is a measure of the local mass-flow capacity of the differential stream 
lamina bounded by the adjacent stream surfaces of constant t) and of constant 77 + 677. 

From equations (LOa) and (10b), for equal values of d^ and 677, 


A 

B 



(10c) 


which is the ratio of sides for the elemental flow area in figure 6. Also, from equations 
(LOa) and (10b), 
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AB = p2j^^sit.20 
dip di} 


or, from the continuity equation (9d), 


From equations (9d) and (10b), 


^B = £_sin_e 

q 


^=qB 

dn 


and equation (9i) becomes 


sin'^e 

Likewise, from equations (9d) and (10a), 


Vip = ■ ” ®3 COS 6^ 

\ / 


-^=qA 

dm 


(lOd) 


aoe) 


ao£) 


(log) 


and equation (9j) becomes 


Vrj = 


= -2±. ^e- - e 


sln^e 


Cg COS 0 


') 


(lOh) 


Equations (lOf) and (lOh) are final forms for the gradients of ip and 77 , respectively. 
The gradient of cp is given by equation (5g) . 

In summary, from equations (5f), (10 e), and (lOg), 


^ =q 

(5f)l 

ds 

^ = qB 

(lOe) ' 

dn 

-^=qA 

aog) 

dm 



(lOi) 


The first of these equations is the result of irrotational fluid motion; the second and 
third result from continuity. 
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Flow field in transformed (;p, t? - space . - The upstream and downstream bound- 
aries of the flow field in transformed -space can be, for convenience, plane sur- 

faces of constant <p . If the value of (p at the upstream boundary is taken as zero, the 
downstream value (p-^ is, from equation (5f), 



where the subscript D refers to the downstream boundary, and q(s) is the prescribed 
distribution of velocity q as a function of distance s along streamlines on the lateral 
surface of the flow field. If the downstream surface is a plane of constant these 
prescribed distributions of q(s) must be such that the integral in equation (11) is the 
same for all these streamlines. 

The boundary shape of the flow field on planes of constant (p in (/?, e/), 77 -space de- 
pends on (1) the specified shape of the upstream boundary surface in x,y,z-space and 
(2) the somewhat-arbitrary, specified shape and orientation of constant ^ and constant 
77 lines on that upstream boundary. However, whatever the shape of the boundary on 
the upstream (^-plane in </?, ^,T7-space, the same shape exists for all planes of constant 
cpy because the paired values of ^ and 77 for each streamline on the lateral boundary 
remain constant. Thus the flow field in <;p, 77 -space is a right cylinder (fig. 9) or a 
rectangular parallelepiped. 

Outline of method . - When q is specified as a function of s along streamlines on 
the lateral surfaces of the flow field, then, from equation (5f) and with uniform veloci- 
ties on the upstream and downstream planes, the velocity q is known everywhere on 
the boundary in (7?, 77-space. For these known boundary values of q, a governing, 
partial -differential equation, derived from the laws of potential flow expressed in cur- 
vilinear coordinates, is developed for the distribution of q within the flow field in 
^,T7-space. And from this distribution of q within the field, auxiliary equations are 
used to compute flow directions and the shape of the flow field in x, y, z-space. 


Governing, Partial- Differential Equation for Distribution of 
Velocity in Transformed <77, ^,77 -Space 

The governing, partial-differential equation for the distribution of velocity q in 
transformed 77-space is developed from the conditions of continuity and irrotation- 
ality for three-dimensional, compressible flow. This general equation is then reduced 
for the special cases of (1) incompressible flow, (2) plane two-dimensional flow, and 
(3) axisymmetric flow. 
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Vector differential operator V. - Expressed in terms of the curvilinear ^,^,77- 
coordinates, the vector differential operator V is given by 

V = V(^ — + V?/) ~ + Vt] — (12a) 

b(p drp drj 

where the vector gradients V<p, Vip, and V17 are given by equations (5g), (lOf), and 
(lOh), respectively. Thus 

V = e-q — +/e - Co cos ^ + /"e^ - e„ cos ^ (12b) 

1 3^. V 2 3 V 3 2 ^ 


Equation (12b) gives the vector differential operator V in terms of (1) the partial dif- 
ferentials with respect to cpt ipt and tj and (2) the unit vectors e^, and e^. 

Irrotationality condition . - For irrotational flow, the vorticity of the velocity vec- 
tor is zero. Thus 


V X q = 0 


asa) 


or, from equations (lb) and (12b), 


e- -^ + - Co cos 0^ — ^ — + /co “ cos 0^ ^ 

1 acp V 2 3 y V 3 2 y 


sin 0 


sin^0 ^ 


X e^q = 0 


so that, with (see fig. 5) 


^2 ” ^3 ^ 

sin 0 


X ei - eg 


(13b) 


63 - Cg cos 0 

sin 0 


^ ®1 ^ ®2 


the vector equation for zero vorticity is 


(13c) 


e x-^ — e ^ A ^ ^ 

1 Q<p sin 0 dip V ^ - 2 


A - a In q 


— — + — e, 
dip sin 0 8 t| 


A /— -- \ 

+ — - — ■“ eg cos dj X — 1- = 0 


sin^0 


(13d) 
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Because the vectors 8 e^/ 9 ^, 9 e’^/ 8 ?^, and 8 ^-j^/ 8 tj are normal to e^, they lie on 
the plane tangent to the <p surface, as do the unit vectors e^ and e^ (see fig. 5). 
Thus the first, second, and fourth terms of equation (13d) are vectors that lie on the 
tangent plane and are normal to the third and fifth terms are vectors that are nor- 
mal to the tangent plane and, thus, parallel to e^^. Therefore the vorticity component 
normal to the 97 surface, which vorticity component is zero, is obtained from the dot 
product of the first, second, and fourth terms with the unit vector e^ ^ ^ 2 ’ vec- 

tor is normal to that surface (see figs. 3 and 5). Thus 



which, from equation (9h) and from page 18 of reference 5, reduces to 


e 


2 


_ B - — = 0 

dcp dip 


(13e) 


Equation (13e) is the first of three scalar equations resulting from the ir rotational ity 
condition. Each term in this equation is equal to l/q times the local geodesic curva- 
ture of the -coordinate line on the 17 stream surface in x, y, z-space (eq. (B7), 
appendix B). 

Likewise the vorticity component normal to the ip surface, which vorticity cQmpo- 
nent is zero, is obtained from the dot product of the first, second, and fourth termis-in 
equation (13d) with the unit vector e^ x "e^, which vector is normal to that surface (see 
figs. 3 and 5). Thus 

X - _S_ e„ . (e„ >< 

d(p sine diP Sin 0 9r] / V V 

which, from equation (9g) and from page 18 of reference 5, reduces to 



= 0 asf) 

877 d(p 

Equation (13f) is the second of the three scalar equations resulting from the irrotation- 
ality condition. Each term in this equation is equal to -l/q times the local geodesic 
curvature of the (;p -coordinate line on the ip stream surface in x,y, z-space 
(eq. (B14), appendix B). 
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Finally, the vorticity component normal to the (^-surface, which vorticity compo- 
nent is zero, is obtained from the dot product of the third and fifth terms in equation 
(13d) with e^j^, which vector is normal to that surface. Thus 


B 

sin^0 



eg cos 6] 




^ Cg cos 01 X 


sin ^0 


ae. 


ar? 


= 0 


which, from equations (9g) and (9h) and from page 17 of reference 5, reduces to 


_ 9e _ 8e- 

Be., * — ~ - Ae„ • — — 

dip ^ dr] 


- 0 


(13g) 


Equation (13g) is the third of the three scalar equations resulting from the irrotation- 
ality condition. Because q is normal to the cp-surface, this equation does not involve 
the velocity magnitude q. To determine the physical significance of this equation, con- 
sider the vector component itself, rather than its scalar value. The third and fifth 
terms of equation (13d), after equations (lOf) and (lOh) are introduced, become 


ae., ae- _ 

X — i + Vrj X — i = Vq X e- =0 (13h) 

dip dr] ^ ^ 

where Vg x e^ is the surface rotation of on a constant (^-surface, here expressed 
in curvilinear 17 -coordinates. (Note that the gradients Vip and Vtj both lie on the 
(^-surface.) Thus, from equation (13h), the surface rotation of the unit vector nor- 
mal to the (^-surface is zero. This condition is true for the unit vector normal to any 
surface (ref. 6 , p. 208). 

Surface rotation of normal unit vectors . - Similar to equation (13h), the surface 
rotation of the unit vector x Cg normal to the 7 ]-surface is equal to zero. Thus 

Vg X (e^ X 63 ) = "" (®1 ° 

where V^ip is the surface gradient of ip on the Tj-surface. (The spatial gradient V<p 
lies on the t] - surface.) From equation (10 e), 

^S^ = e ^ (14b) 

^ ^ dn ^ 

Thus, from equations (5g), (14a), and (14b), 
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die X e^) _ a(e- x e„) 

e X ^ + Be„ X — ^ = 0 

^ d(p ^ dip 

Expanding the triple vector products and taking the dot product of the result in the di- 
rection of e^ X yield 



or introducing equation (9h) gives 

9e„ ^ \ 9e, 

e„ • — - = Bf cos • — - (14c) 

d<p ^ /dip 

In the same fashion, because the surface rotation of the unit vector e^ x e^^ normal 
to the i/)-surface is zero. 


— /— • — \ 

e„ * Af e - e„ cos 0)* (14d) 

dcp \ ^ /dr) 


Equations (14c) and (14d), which result from zero surface rotation of the unit vec- 
tors e- X e_ and x e- normal to the 17- and jp-surfaces, respectively, can be com- 

JL /J O X ^9 

bined as follows. From equation (8b) and the fact that sin 0 + cos‘^0 = 1, 


— ^^3 , — ^^2 _ sin^0 8 In sin 9 

e , -h 

dcp ^ dcp cos 9 dcp 
or, from equations (14c) and (14d), 


B 

sin^9 



e^ cos 





dip 


A 

siifie 



\ 

e„ cos 0) • — i- = 
^ / dr) 


Be, 


8 e. 


dip 


+ Ae. 


^^1 ^ 8 In sin 9 
dr) d(p 


Equation (L4e) is used later to simplify the continuity equation. 


a4e) 
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Continuity condition . - For steady, compressible flow, the divergence of the product 
of the velocity vector and the fluid density is zero. Thus 


V ♦ (pq) = 0 (15a) 

or, from equations (lb) and (12b), 



which expands to give 


8tn , q ^ alnp B /- 
9'?’ sln^fl ^ ^ 


COS 01 




-A “ 

6o - COS 0 
9 V O Z 

sin^0 


\ 

9p = 0 

/ dn 


Introducing equation (14e) gives the continuity equation in transformed (p, ?/), 77 -coor- 
dinates; 


9is_a + lbLa+B?„ 

dtp b(p ^ 


8e- _ 9e., 9 In sin B 

— + Ae„ * — — 0 

dip dr] dtp 


(15b) 


Surface divergence of flow in stream laminae . - In a steady, three-dimensional, 
compressible flow field, consider the fluid in a differential stream lamina bounded by 
adjacent stream surfaces of constant 77 and of constant 17 + dry. In general, the ele- 
mental thickness dm sin 6 (fig. 8) of this lamina varies from point to point on the 
77-surface of the lamina. From continuity, the surface divergence of the flow 
pq dm sin B in the lamina is zero , Thus 


Vg * (pq dm sin 0) = 0 

where, as in equation (14a), 

® dcp ^ dip 

so that, from equations (5g), (14b), and (16b), equation (16a) becomes 


e — + K b • (qe- B drj) = 0 

,^dcp QipJ ^ 


(16a) 


(16b) 


which, with dry constant for the lamina, reduces to 
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(16c) 


Be - 9 In q 3 In B 

^ dip dtp d<p 

Likewise, for the flow in a differential stream lamina bounded by adjacent stream 
surfaces of constant ip and of constant ip + dip, 

Ae„.^= (16d) 

dr] d(p d(p 

Equations (16c) and (16d) are the continuity equations for flow in the two families of 
stream laminae. Each side of equation (16c) is equal to 1/q times the local geodesic 
curvature of the ^-coordinate line on the rj stream surface in x, y, z-space 
(eq. (B8), appendix B); likewise each side of equation (16d) is equal to -l/q times the 
local geodesic curvature of the rj-coordinate line on the ip stream surface 
(eq. (B13), appendix B). 

If equations (16c) and (16d) are entered into the general continuity equation (15b), 

9 In p 9 In q 9 In A _ 9 In B ^ 9 In sin 0 _ ^ 
dcp dcp dcp d(f dcp 


or 


P ^ = constant 
qAB 

where, from equation (lOd), the constant is 1. 

Governing differential equation . - The governing differential equation for the distri- 
bution of velocity q in (p,ip,r]-spa.ce is derived from the continuity equation (15b) and 
from the two equations for the irrotationality condition, equations (13e) and (13f). 

As already observed, any three-dimensional, steady-flow field is defined by the 
spatial distribution of the velocity q and its direction where, from equation (la) 
and/or figure 1, the unit vector e"^ is defined by the three angles ci^, and 
From equation (la), the three angles are related by 

— _ 2 2 2 

e^ • Cj^ = 1 = cos + cos + cos y^ (17a) 

Thus, when equation (17a) is taken into account, the flow field is defined by three un- 
knowns; the velocity q and two flow angles. The solution, therefore, involves three 
equations . These equations are the continuity equation (15b) and the two equations for 
the irrotationality condition, equations (13e) and (13f). All three of these equations 
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involve partial derivatives of q and only. After appropriate second partial deriv 
atives are obtained, the three equations are combined to eliminate the second partial 
derivatives of e^^. 

Thus, from the continuity equation (15b), after the partial derivative with respect 
to (p is taken, 


In q ^ 9^ In p ^ 
d(p^ dcp^ 


d In sing 
dcp^ 


+ B 


d<p 


, . 9 In A — 
+ A e. 

dcp 



9e, 

9e„ 

9e_ 


. l + B 

2 

^ 1 

2 

dip 

dcp 

dip 

9e. 

9e„ 

9e_ 


1 

+ A — ^ . 

1 

+ Ae„ 

dr] 

dcp 

dr] 

o 


+ Be, 


9^e^ 
dcp dij) 


9%. 


dcp dr] 


= 0 


(17b) 


where, from the irrotationality equation (13e), after the partial derivative with_respect 
to ip is taken, 


;e . = _ B ^ ^ + b 2 ^ ^ 3 In q g2 3^ In q 

^ dcp dip dip dcp dll 34 g . 2 


(17C) 


and where, from the irrotationality equation (13f), after the partial derivative with re- 
spect to T] is taken. 


= _ A ^ ^ 9 1"^ 9 Inq ^ ^2 3^ In q 

^ dcp 8t) dri dcp 3t) 9t) g 2 


(17d) 


Therefore, from equations (17b), (17c), and (17d), after equations (16c) and (16d) are 
introduced, the governing differential equation for In q becomes 


dcp"^ dcp‘ 


dcp‘ 


8^ In q . 8^ In 0 , 8^ In sing . ) , _y 

y9(p dip dip dcp^ 


9eg 9e^ 9e^ 9e^ 

dcp 9tj drj dcp 


-i- I Q ^ 9 In q ^ 9^ In q \ + ^2 / 9 In A 9 In q ^ 9^ In q ' 


dip dip 


dip' 


drj dr] 


dr]' 


9 In B [ 9 In q ^ 9 In B \ 9 In A / 9 In q ^ 9 In A 

dcp \ dcp ^(p / 9^ V 9^ 9(^ 


)=0 


a7e) 
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Equation (17 e) can be further simplified by noting (appendix C) that the fourth and fifth 
terms are related to the 
faces, respectively, by 


terms are related to the total curvatures and of the 77 and ^ stream sur- 


2 ^ 9^ 


a7f) 


K 

( 

q' 


9e„ 0e, 8e„ Be. 


,2 \ d(p dr} drj 


dcp 


(I7g) 


Thus, from equations (17 e), (17 f), and (17 g). 


B^ In q , B^ In p , 9^ In sin 6 ^77 


dcp‘ 


dcp* 


dcp‘ 


2 2 

q q 


+ ( 9 In B 8 In q ^ 8^ In q \ _ 8 In B / B In q , B In B 


dip dip 


dip' 


dip \ dip 


dip 


^2 /a In A B In q , B^ In q ] B In ^ / B In q , 8 In _ 


dr] 


drj 


dr]‘ 


dip 


dtp 


dip 


(17h) 


Equation (17h) is the governing, second-order, partial-differential equation for the dis- 
tribution of In q in transformed ip,ip,r} -space. This equation is also developed (appen- 
dix D) directly from separate considerations of the two-dimensional continuity and irro- 
tationality flow conditions in 77 and ip stream laminae; the results are then interre- 
lated by the simple continuity equation (lOd) for three-dimensional flow. 

Incompressible flow . - For incompressible flow, the density p is 1.0. Thus the 
second term in the governing differential equation (L7h) disappears, and the density no 
longer affects the values of A and B given by equations (10a) and (10b), respectively. 

Plane two-dimensional flow . - For plane, two-dimensional flow. 


V 




0 = 90' 
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so that, from equations (10a) and (10b), 


A = p^ sin 0=1 


d?/; 

dm 

d?7 


B = p sin 6 = p 


from which 


8 In A / 8 In q 8 In A \ _ _ 


dcp 


d<p 


dcp 



8 In q _ 8 In q \ _ ^ 
dcp dcp / 


Thus, for plane, two-dimensional flow, the governing differential equation (17h) reduces 
to 


8^ In q , 3^ In p , ^2 / 3 In p 3 In q , 9^ In q\ 3 In p /d In q , 8 In p\ _ q 

a<p^ dq}^ \ ^ ^ 94,^ J Sq> \ a<p a<p / 


Equation (18a) is the governing differential equation (7) in reference 1. 
- For incompressible flow (p = 1), equation (18a) becomes 


3^ In q , 3^ In q _ g 
d<p^ dip^ 


(18b) 


Axisymmetric two-dimensional flow . - For axisymmetric, two-dimensional flow, 

K =0 ^ 

-n 


6 = 90® 


r 


(19a) 


(continued) 



dm _ ^ 
dr) 


and, except for 9eo/9n> 


dn _ 1 
6ip prq 


-9 =0 


y 


dn 


J 


where r is the radius from the axis. Thus, from equation (LOa), 


A dn r\ 3- 
A = p — sin 0 = — 


d^ 


rq 


and, from equation (10b), 


8 In A ^ _ /a In q 8 In r 

8(^ \ dcp d(p j 


8^ In A _ In q , 8^ In r 
d(p^ \ d(p^ d(p' 


B = P sin 0 = pr 
dTi 


8 In B ^ 8 In p ^ 8 In r 

d(p dcp dcp 

8 In B ^ 8 In p 8 In r 

dip dip dip 


Also, from equations (C6) and (19b), 








(19a) 


(19b) 


(19c) 


9^ In r 8 In r j 8 In q ^ 8 In r^ 


2 ^2 
q dcp 


(19d) 


dcp \ dcp 9(^ y 

Thus, from equations (19a) to (I9d), the governing differential equation (17 h) reduces to 
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3^ In q , 9^ In p , 9^ In r 
d(p^ dcp^ d(p^ 


liiLa + 

dcp 


9 In r 
dcp 



8 1np I 

dcp 


9 In r 
dcp 


^22 
+ p r 


( 9 In p ^ 9 In r\ 

dip dip / 


il!La + 

dip 


d^ In q 


dip^ J 


= 0 


(I9e) 


For axisymmetric, two-dimensional, compressible flow, equation (19e) is the govern- 
ing differential equation for the distribution of In q on the transformed cp^ j,ii-plane. 

For incompressible flow (o = 1), equation (19e) reduces to 


9^ In g 
dcp^ 


+ 9 In r _ 9 In r [ 9 In q ^ 
2 dcp \ dcp 


9 In r ^ ^ ^2 [ 9 In r 9 In q 


dcp^ 


dcp 


dip dip 



(19f) 


Equation (19f) is the governing differential equation (11) in reference 2. 


Construction of Physical Flow Field 

For prescribed distributions of velocity along the boundary streamlines, the ob- 
jective of this design method is the construction of the physical flow field, particularly 
the shape of its boundary, in x,y,z-space. From every finite-difference solution of 
the governing differential equation (17h), in a series of iterative solutions, each with 
fixed distributions of A, B, p, and sin 6 throughout the flow field, the distribution of 
velocity, that is, of In q, is known in <p, i 7-space. From this distribution of velocity, 
the shape of the flow field in physical x,y, z~space can be determined by methods, de- 
veloped in this section, using various equations introduced during the development of 
the governing equation (17h). 

Outline of construction method . - After the governing equation (17h) has been 
solved for the distribution of In q (by finite-difference techniques), the unit vectors e^, 
'62, and e^g; the physical coordinate x, y, and z; the angle 0; and the parameters A 
and B can be determined as follows: 

(1) Assume that e^^, e2» x, y, z, 0, A, and B are known for a given ^-sur- 

face. For example, this surface is initially the upstream boundary (cp =0.0), where, 
as shown later, x, y, and z are known and with q equal to e^ (q = q^ =1.0) and 0 
equal to 90®, 
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A = 1 
B = 1 








( 20 ) 


(2) The object then is to determine e^, 62, e^, x, y, z, 0, A, and B for the next 
downstream (^-surface. 

(3) This object is achieved by the following procedure, which is described more 
fully in subsequent sections: 

(a) From equations (7a), (13e), and (13f), determine the distribution of on 
the next potential surface in x,y, z-space. 

(b) From equations (7c), (13e) (with e^ • 62 = 0)» and (14c), determine the dis- 
tribution of 62 on the next potential surface . 

(c) From equations (7e), (13f) (with e^ ■ e^ = 0), and (14d), determine the dis- 
tribution of on the next potential surface . 

(d) Determine the distribution of 9 on the next potential surface from equa- 
tion (8 c). 

(e) Determine the distributions of x-, y-, and z-coordinates on the next poten- 
tial surface by methods to be presented. 

(i) Determine the distributions of A and B on the next potential surface from 
equations (10a) and (10b), respectively, in which p and sin 6 are known and where 
6n/6ij) and dm/drj follow directly from the distributions of x, y, and z obtained in 
step (e). (Alternative methods, which are discussed later, may be used for the deter- 
mination of A and B . ) 

Unit vector e^. - From equation (7a), 

— — 2 2 2 

• e^ = cos + cos /3^ + cos 7^ = 1 (21a) 


8e, — 8 cos q;, — 8 cos — 8 cos y, 

— i=i L+j i+k U (21b) 

d<p d(p dcp dcp 

where, from equations (13e), (13f), and (2 la). 
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in which 


cos cos ( 3 ^ 0 cos ( 3 ^ COSYj^ 0cosy^ 


d(p 


COS a. d(p cos dcp 


(21C) 


0 COS (3^ ^2^3 


d(p 


^1^4 - ^2^3 


(21d) 


* * 

0cosy^ ^2^1 ^ 


1 '^V2 


d(p 


P P - P P 


(21 e) 


* = B 8 i”q 


C^ = B 


dip 


% 

C2=^ A 


=. 

dr] 


>> 


(2lf) 


P = cos q;„ cos jS - cos a, cos /3r 


?2 = cos cos - cos o;^ cos /3^ 


Pg = cos q ?2 7i~ 72 


P^ = cos Of 2 cos 7 -^ - cos cos y ^ 


(2lg) 




Integration of equations (21c), (2ld), and (2le) determines the distributions of cos a^t 
cos /3^, and cosy^, respectively, and, thereby, from equation (7a), the distribution 
of e^ on the next potential surface. 

Unit vector e2- -- From equation (7c), 


— 2 2 2 
Cg = cos Qfg + cos ^2 + cos y2 = 1 


(22 a) 


06o _ 0 COS Q!^ _ 0 COS /3o _ 0 COS y o 

—2 = 1 l+j 2^,^ r _2 (22b) 

dcp dcp dcp dcp 

where, from equations (13e) (with (I'^c), and (22a), 
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(22c) 


acosQjg cos /?2 9 cos /?2 cos 72 9 cos 72 

d(p cos ttg 9 ^ cos Qfg 9 ^ 

8 cos “2 ^ ^3^8 ^2 

^ ^ 6^7 " ^ 5^8 

8 cos 72 ^2 ^3^6 ^2 

^ 5^8 ■ ^ 6^7 

in which 


(22d) 


(22e) 


* ^ 
S = ® 


8 cos O'- 8 cos / 3 - 

(cos " cos o?2 cos B) ^ + (cos jS^ - cos ^2 — 


L 




+ (cos 7 2 - cos 72 cos 6) 

Pg - cos cos ^2 ■" ^2 ^3 

Pg = -P^ = cos cos ^2 ~ ^2 ^1 

P^ - cos Og cos 72 “ cos Q'2 cos 7^ 


pg = -Pg = cos o?^ COS 72 - cos Q'2 cos7j^ 


dip 

d cos 7 ^ 
dip 


(22f) 


(22g) 


Integration of equations ( 22 c), ( 22 d), and ( 22 e) determines the distributions of cos 0:2, 
cos ^2* cos 721 respectively, and, thereby, from equation ( 7 c), the distribution of 
■62 on the next potential surface. 

Unit vector e^. - From equation ( 7 e), 

_ 2 2 2 

03 • eg = cos Q!g + cos jSg + cos 73=1 ( 23 a) 


8e„ 8 cos Ckfo -.9 cos _ 9 cos 

i + j i+k ^ (23b) 

d(p d<p d(p d(p 
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where, from equations (13f) (with e^ * e^ = 0), (14d), and (23a), 


9 cos Q!g cos /3g 9 cos cosy^ 90087 ^ 


d(p 


cos a. 


d(p 


cos a. 


dcp 


9 cos jSg ^2^9 ^3 ^4^10 ^3 


d(p 


p p _ p p 

12^9 ^ 11-^10 


9 cosyg ^2^11 ®3 ^4^12 ®3 


9</? 


p p _ p p 

11 10 12-^9 


(23c) 


(23d) 


(23e) 


in which 


* 

C^ = A 


9 cos O' 9 cos /3, 

(cos ^2 - cos cos 6) + (cos ^2 " 


9n 


9?i 


9 cos y 

+ (cos y - cos y „ cos 0 ) 

^ ^ J 


(23f) 


Pg = -P^ = cos 02 cos 73 “ cos Og cos 72 
P^Q = -P^ = cos cos 7 2 - cos o^ cos y^ 




(23g) 


Pll = -Pg = cos 02 Ps ” ^3 ^2 

^12 ” ”^2 ~ ^1 ^3 " ^3 ^1 _> 

Integration of equations (23c), (23d), and (23e) determines the distributions of cos o^, 
cos jSg* and cosy^, respectively, and, thereby, from equation ( 7 e), the distribution of 
e^ on the next potential surface. 

Parameter A . - The parameter A comes directly from its definition, equa- 
tion (10 a): 


A 


= p sin 9 

Aip 


(24) 
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where An follows directly from known values of Ax, Ay, and Az (corresponding to A^ 
and obtained as described in a later section) on the next potential surface. Equation ^4) 
determines the distribution of parameter A on the next potential surface. 

Parameter B . - The parameter B comes directly from its definition, equa- 
tion (10b): 

B = p — sin e (25) 

At? 

where Am follows directly from known values of Ax, Ay, and Az (corresponding to 
At? and obtained as described in the next section) on the next potential surface. Equa- 
tion (25) determines the distribution of parameter B on the next potential surface, 

Although in theory the values of A and B are determined by equations (24) and (25), 
respectively, in practice, the values may be better determined, during the iteration 
process, by equation (lOd) and the ratio of equations (24) and (25) or by equation (lOd) 
and the difference of equations (16c) and (I6d), which difference gives 



Coordinates of flow field in x, y,z-space . - The x, y, z-coordinates corresponding 
to points in the flow field of (p, 2 p,r?-space, which coordinates establish the shape of the 
flow field in physical x, y, z-space, are obtained from the direction cosines. Thus, 
from equation (7b), along a ^-curve of constant ip and tj. 


9x 

0s 


cos 


dx d(p 
dcp ds 


so that, from equation (5f), 

/ (p 

/cos a A 

„ [-7-r 

Likewise 

dy _ dy ^ ^ 

— ^ = -^ q = cos (3^ 

ds dcp ds dcp ^ 


(26a) 


so that 
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(26b) 


And likewise 


so that 


8z 9z 9z ^ 

9s d(p ds d(p ■*■ 



(26c) 


Equations (26a), (26b), and (26c), which are integrals along streamlines of constant ip 
and T], are sufficient to determine the physical shape of the flow field in x,y, z-space. 
In particular, for streamlines along the boundary, these equations determine the shape 
of the boundary . 

In addition, along ip- and r]-curves on any potential surface of constant cp, the 
X, y, z-coordinates can be determined by other, redundant integrals. Thus, along a 
ip-curve^ 


9x 

9n 


= cos a. 


dx dip 
dip dn 


cos 0!2 


so that, from equation (lOe), 


X 



(26d) 


where the subscript x refers to conditions at the lower bound of the integral. Like' 
wise 


9n dip dn dip ^ 


so that 
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y = y^ 


(26e) 


And likewise 


so that 



^ = ^ a = ^qB = coS7, 

9n dip dn dip 



(26f) 


Equations (26d), (26e), and (26f) determine the x,y, z-coordinates on any potential sur- 
face of constant <p from integrals along lines of constant rj on that surface in cp yip, r]- 
space. 

Also, along an T/-curve on any potential surface, 


9x 

9m 


= cos Of. 


a„ 

dr] dm 

so that, from equation (lOg), 



Likewise 

= COS /3o 

dm drj dm dr] 


(26g) 


so that 
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( 26 h) 



And likewise 


y = y^ 



dT7 


8 z 9z dn 9z * 

^ = — qA - cos y „ 

9m dr) dm drj 


so that 



(26i) 


Equations (26g), (26h), and (26i) determine the x, y, z-coordinates on any potential sur- 
face of constant q> from integrals along lines of constant on that surface in (p,iptri- 
space. 

The shape of the flow field in x,y, z-space, including the boundary shape, is given 
by the x, y, z-coordinates determined from equations (26a) to (26i). 


DESIGN PROCEDURE 

The design procedure, using equations developed in the section THEORY OF DE- 
SIGN METHOD, is outlined in this section. After general aspects are considered, spe- 
cial problems associated with internal flow fields and external flow fields are discussed, 
and a brief step-by-step numerical procedure is presented. 


General Considerations 

General aspects of the design procedure that are independent of the internal or ex- 
ternal nature of the flow field are considered here. These aspects include (1) the flow 
field configuration in (^, -space, ( 2 ) finite-difference grid points of the flow field in 
cpfip,r)-spacey (3) the finite difference form of the governing differential equation, (4) the 
specified velocity distribution on boundaries of the flow field, (5) stagnation points, and 
( 6 ) planes of symmetry. 

Flow field configuration in 77 -space . - Because the lateral boundary of the flow 
field in x,y, z- and in 17 -space always contains the same streamlines of constant ^ 
and 17 , the shape of every potential plane in 77 -space is the same, and because the 
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upstream and downstream boundaries are surfaces of constant (p, the flow field config- 
uration in transformed (^, jp, r)- space is a right cylinder, as shown in figure 9, or a rec- 
tangular parallelepiped. 

The shape of the flow field on the potential planes in (p, ?/>, 77 -space depends on the 
specified shape of the upstream boundary area in physical x, y, z-space and on the spe- 
cified arrangement of the fluid lines on that area (lines of constant ip and t), see the 
section Stream functions ^ (x, y , z) and t] (x, y, z)) . Because this arrangement is some- 
what arbitrary provided only that equation ( 6 d) is satisfied), it is generally convenient 
to have the fluid lines of constant ^ intersect the fluid lines of constant 77 at right 
angles. Thus the angle 0 is 90^ at the upstream boundary. As a result, from the 
continuity equation (9d), at the upstream boundary, where p and q are 1.0, 

dip dr] - dn dm (27 a) 

In addition, if a Cartesian grid of fluid lines is used (as opposed to a polar grid, for 
which 6 is also 90^), then, because the flow is uniform everywhere on the upstream 
boundary. 


_ d77 
dn dm 

From equations (27 a) and (27b), it follows that 

dip = dn 
dr 7 = dm 


(27b) 


(27c) 


Thus the shape of the upstream boundary in cp ,ip,r}~ space is the same as the specified 
shape of the upstream boundary in x,y, z-space. 

If a polar grid of fluid lines, instead of the Cartesian grid, is used at the upstream 
boundary in x,y, z-space, the shape of the upstream boundary in <;p, 77 -space is differ- 

ent. This case is considered later in the report. 

Grid points of flow field in cp , ip , r]-spSLCe . - The finite-difference solution of the 
governing differential equation (17h) requires a three-dimensional grid of finitely 
spaced points in transformed (p,ipfT] -space; at these grid points. In q and other proper- 
ties of the flow field are determined. This grid is formed by the intersections of planes 
of constant (p, ip^ and 77 spaced a^, a 2 » and a^ distances, apart, respectively, as 
shown in figure 9. These distances, although constant between any two adjacent planes, 
need not be the same for all planes. The grid also includes the intersections of ip 
and 77 stream planes with the boundary of the flow field on potential planes of con- 
stant cp (fig. 9) . It is at these boundary grid points that the velocity distribution is 
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specified. Obviously the Ug or spacing of a boundary grid point from its adjacent, 
internal grid point on the potential plane is not, in general, uniform. 

From the above discussion, for every internal grid point of the flow field in (^,zp, 97- 
space, there exists a three-dimensional, finite-difference star (fig. 10) consisting of 
the grid point itself (marked 0 ) and its six adjacent grid points (marked 1 to 6). In this 
star, the lengths a^, a2» . . are not, in general, equal. It is this star that is 

used in the finite-difference solution of equation ( 17 h) . The star is in its simplest 
form, and, in some applications, greater complexity, involving more grid points, may 
be required. 

Corresponding to the grid in 77-space is a grid in x,y, z-space. This grid in 
x,y,z-space constitutes the solution and is given by x, y, z-coordinates determined, af- 
ter the solution of equation (L 7 h), by equations ( 26 a) to ( 26 i) at the grid points in <p,2p-rf- 
space. 

Finite-difference form of governing differential equation . - In the governing differ- 
ential equation ( 17 h), first and second partial derivatives of in q with respect to cp, if) ^ 
and 17 can be expressed in finite-difference form by assuming that, in the immediate 
neighborhood of the central grid point 0 in the finite-difference star of figure 10, the 
distribution of In q can be represented by a Taylor series expansion, which, for the 
three-point systems of this star, is equivalent to a parabolic distribution of In q 
through the three grid points in each of the three directions (p, ip, and 17. For such 
distributions, the first derivative at the central point 0 is given by 



Q+ ~ Qq 


Ax, 

(Ax + Ax_) 

Ax 


Q_ ” Qq 

Ax 

(Ax^ + Ax J 


(28 a) 


where Q is In q, x. is distance in the cp-y ip-, or 17 -direction, the subscripts + and - 
refer to values in the positive or negative direction from the central grid point 0, and 
Ax is the distance between adjacent grid points in cp,ip,rj-spa.ce. The second partial 
derivative is likewise given by 


- Qq) 2 (Q_ - Qq) 

Ax^ (Ax^ + Ax) Ax_ (Ax^ + Ax_) 

If the distances between adjacent grid points in a given direction are equal, 

Ax_j_ = Ax_ = Ax 



( 28 b) 
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and equations (28a) and (28b) reduce to 



_ Q+ + Q- ~ ^Qp 

l^axf / ^x)^ 

\ 1 /q 


(28c) 


(28d) 


Introducing the finite-difference equations (28a) and (28b) into the governing differ- 
ential equation (17h) and using the nomenclature from the finite-difference star in fig- 
ure 10 result in 






CqQo iR 


(28 e) 


where 


_ 9^ In p I 9^ In sin 6 


dcp^ 


d<p* 


_ /g In bV _ /d In aV _ s 

\ d(p J y dcp j exp (2 Qq) 


^1 = 




9 In B 9 In A 


d(p 


dcp 


C = ^5^^ ^ 9 In b\ 

^ “2<^2 + \^5 ^ / 

c - ^6^ /_2_ 9 In a\ 

^ “3<^3 * V W ^ ) 

; / 2 , 9 In B 9 In A^ 

^ ^4(^1 + d<p d(p 




y (28f) 


2 

c = ^2^ /_2 _ 9 In b\ 

“ ='5<“2 + V ^ / 

c ^3^^ _ 9 In A\ 

Eg (Eg + Eg) 9rj y 






+ B‘ 


9^ 


‘2*^5 



1 + A' 


2/9 I n A ^6~ ^3 , 2 


9rj 


^3^6 



^3^6; 
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and (ft is the residual error, which measures the degree to which the values of Q at 
the grid points of the finite-difference star in figure 10 satisfy the governing differential 
equation. The value of (ft is reduced, or eliminated, at each interior point of the 
three-dimensional grid in ^,?/>,T 7 -space by iteratively changing the values of while 

maintaining constant values for the coefficients given by equation (28f) . After the resi- 
duals at every interior grid point have been reduced to sufficiently low values, new val- 
ues for the coefficients are computed from equations (28f), and the iterative process is 
repeated. This procedure is continued until the values of the maximum residuals have 
converged within the derived accuracy; at which point, the governing differential equa- 
tion (17h) has been solved for the specified distribution of Q ^.e. , In q) on the bound- 
ary of the flow field in 77 - space. 

Specified velocity distribution on boundary of flow field in 77 -space. - As de- 

scribed in the section Preliminary Considerations, the specified distribution of velocity 
q as a function of path length s along streamlines on the laterial boundary of the flow 
field in x, y, z-space translates directly into the distribution of q as a function of veloc- 
ity potential cp along streamlines on the lateral boundary of the flow field in 77 - 
space. This direct relation results from equation (5f). hi this section, therefore, con- 
siderations are limited to the velocity distribution as a function of cp on the boundary 
of the flow field in (^, 7j-space. 

For proper use of the governing differential equation (17h), the upstream boundary 
(cp = 0 ) and the downstream boundary {<p = cp-^) should be sufficiently far upstream and 
downstream, respectively, from the region of the lateral boundary in which velocity 
changes are specified that the velocity is constant, as specified,, not only on both the 
upstream boundary (where q - 1 . 0 ) and the downstream boundary (where q = qj))» but 
also over substantial regions of the lateral boundary adjacent to the end boundaries. 

In the region of the lateral boundary where velocity changes are specified, these 
changes should be at least piecewise continuous, both in the direction of cp along the 
streamlines and in the direction tangent to the boundary and normal to the streamlines . 
To assure adequate smoothness in the variation of prescribed velocity on the laterial 
boundary, specify the velocity q as a function of cp and as a function of the parameter 
P|j, where p^ is path length along the prescribed perimeter of the upstream flow area. 
Thus, on the lateral boundary. 


from which, because 


q = 

(29 a) 


(29b) 

'n ='^(Pu) 

(29c) 
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it follows that 


q = q(<p,2/>,77) (29d) 

Stagnation points . - Stagnation points, at which the velocity q is zero so that In q 
becomes negative infinity, are singularities that occur on the boundaries of all bodies in 
external flow fields and, for internal flow fields, at the junction of branched ducts or on 
the nose and tail surfaces of centerbodies within ducts . The mathematical problems 
associated with these singularities are not treated in the present report . Instead low 
but reasonably positive velocities can be specified on the boundary in regions that would 
otherwise include a stagnation point . As a result, these regions of the boundary are 
cusped in physical x, y, z-space. 

Planes of symmetry . - Planes of symmetry, for which the flow fields on each side 
are mirror images, can occur in both internal and external flow fields. These planes 
of symmetry are either parallel to or normal to the velocity vectors on the plane . Such 
planes can be treated as boundaries of the flow field in both cpyip,7]- and x,y, z-space. 

On these boundaries of planar symmetry, the velocities are not specified, but instead 
the derivative of the velocity normal to the boundary is zero. 


Internal Flow Fields 

Internal flow fields are ducted flow fields in physical x, y, z-space. Examples are 
simple ducts, branched ducts, annular ducts, and ducts with centerbodies, including the 
nose and/or tail region. The three-dimensionality of these ducts implies that they are 
unsymmetrical, although the ducts can have planar symmetry (q.v.). 

Upstream boundary . - The upstream boundary is a flat surface of constant velocity 
potential (<^^ = 0,0) in both x, y, z- and cp, 17-space. The orientation of this boundary 
in x,y, z-space is arbitrary, but in this report, it is taken as the y, z-plane of constant 
X equal to 0.0, as shown in figure 11. 

The shape of the upstream boundary in x, y, z-space is an arbitrary input that is 
specified by 




(30 a) 


Zu “ (30b) 

where p^ is path length along the perimeter starting from the point marked on the pos- 
itive y-axis in figure 11. The area of the upstream boundary, which area is equal to 
1.0 (see the section Nondimensional forms) , is then given by 
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1.0 =i ^ (y^dZy- Zjjdyy) 




(30c) 


If the upstream boundary includes an inner perimeter, such as would be the case for an 
annular duct, the line integral includes the inner perimeter, which is connected to the 
outer perimeter by a line cut, as shown in sketch (a) , 



The intersections of the two sets of stream surfaces, z^(x,y, z) and t](x, y, z), with 
the upstream boundary <^(0,y,z) = 0 in x,y,z -space are two sets of fluid lines that can 
be viewed as generating the stream surfaces as they move downstream with the flow. 
The specified (Input) arrangement of these sets of fluid lines, that is, sets of constant 
^ (0»yu» sets of constant ? 7 ( 0 ,y^, z^)* on the upstream boundary is somewhat 

arbitrary, provided only that equation (6d) is satisfied, in which case the local angle 0 -^ 
with which the line of one set intersects that of the other is ever}nvhere greater than 
zero and less than 180^. However, as observed previously, the arrangement of these 
fluid lines, together with the shape of the upstream boundary, both of which are input by 
the designer, dictates the shape of the flow field in 77 -space. It is, therefore, 
convenient to arrange the fluid lines so that is ever 5 rwhere equal to 90® on the up- 
stream boundary. In particular, if the fluid lines are arranged so that 
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(30d) 


^ = Cy^ 

V “ 

where the constant C is the same for both equations, then from the continuity equa- 
tion (9d), 





sin 9 




or, because p^=1.0, q^ = 1.0. 


and 0 y = 90 


o 


C = 1.0 


from which equations (30d) become 


<^ = yu ^ 


(30 e) 

Equation (SOe) is similar to equation (27c). 

Thus, in summary, the desired shape of the upstream boundary 
specified, by the designer, on the y^, z^- plane, by 

in X, y, z-space is 

1! 


(30a) 

Zu - Z|j(Pu) 


(30b) 

with fluid lines arranged on the upstream boundary so that 


* = yu \ 

t? = zu I 

> 

(30 e) 


where the nondimensional values of y^j and are such that the upstream area 
is 1.0. (In practice, the upstream shape is first specified in dimensional units from 
which the dimensional value of the upstream area is determined; the nondimensional co- 
ordinates, for which the upstream area is 1.0, are then obtained from equations (4a).) 
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For the previously specified arrangement of fluid lines, the arbitrary shape of the 
upstream boundary with a grid of j/j and 17 fluid lines is that shown in figure 11 . It is 
noted that the Aif) grid spacing need not equal the A17 spacing, and particularly, in 
this general case, the spacing for grid points adjacent to the boundary is variable. Be- 
cause jp = y^ and 77 = z^, the grid spacing and the shape of the upstream boundary 
are unchanged by the transformation from x, y, z-space to j/), 77-space; in addition, this 
shape and grid remain unchanged with (p in j/j, 77-space. Of course, for upstream 
boundaries that are square or rectangular, the rectangular grid in figure 11 is ideal, 
because the grid spacing adjacent to the boundary need not be variable. (In this case, 
the grid points at the four corners of the boundary do not appear in the finite-difference 
star (fig. 10 ). Thus, in effect, these upstream boundaries may or may not have rounded 
corners with radii equal to or less than the Ay^j or Az^ grid spacings, whichever is 
smaller. ) 

For circular upstream boundaries in x,y, z-space, including circular annuli and 
sectors of circular annuli, a circumferential and radial arrangement of the ip and 77 
fluid lines, respectively, is more convenient. For this arrangement, »u is again 90^, 
and 


where 


dnu=dru 


dmu = rude^ 




J 


( 30 f) 




0 ^ = tan 


-1 


yu 






( 30 g) 


From equation (30 c), with 


yu = cos 9 \j 




( 30 h) 


the nondimensional area of the upstream boundary becomes 
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(30i) 


or, for the special case of constant values of r^ with varying 0jj, 

l-0~fu,max-4.min) <30j) 

where r^j „ and Ttj are the outer and inner radii, respectively, and A0U is 
U, max U,mm » r- y 

the included angle. 

If, by definition, for this special case of constant r^ with varying 0^, 


U 


A0 


(30k) 


U 


so that the full range of t) values goes from 0.0 to 1.0, then, from equations (30f) and 
from the continuity equation (9d), 


Pxfiv = 


_ 1 /dA / dT) 

\de\ 


or, because p^=1.0, q^ = 1.0, and 0^^ is 90 , 

dip = A0jj r^ dr^j 


U \^VU /U 
o 


from which 


f f 


U 




U, min 


or 


, (2 2 

ip = (rTT - r 


U U, min 


From equation (301), with r^ equal to r^ 


max 


I , f 2 

^max 2 \U,max"^U 


, min^ 


( 301 ) 
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or, from equation (30 j). 


lb =1.0 
^max 

so that the full range of ip values goes from 0.0 to 1.0, as is the case for the 17 val- 
ues . 

For the previously specified arrangement of fluid lines, the shape of the upstream 
boundary is that shown in figure 12, in which the scale of the y^, -coordinates is such 
that the upstream boundary area is 1.0. The fluid lines of constant ip are here 
spaced, according to equation (301), to give equal increments of ip. However, as is the 
case for rectangular Cartesian grids, the Aip and At) grid spacings need not be equal, 
except of course, between adjacent stream surfaces. In any event, for upstream bound- 
aries with combinations of straight- radial and circular- arc contours, the grid spacing 
for points adjacent to the contour can be conveniently constant. Also, as for the rec- 
tangular Cartesian grid (q.v.), because the grid p>oints at the four corners of the bound- 
ary do not appear in the finite-difference star (fig. 10), the upstream flow area can have, 
in effect, rounded corners, unless A0|j equals 360^. However, unlike the Cartesian 
grid arrangement, for which the shape of the upstream boundary remains unchanged by 
the transformation, for this polar grid arrangement of fluid lines, the shap>e of the up- 
stream boundary changes from one with circular features in x,y, z- space to one with a 
rectangular shape in (p, z/), 17 -space. 

Because, for potential flow, the solution is reversible, in some cases the shape of 
the upstream boundary, which shape is arbitrarily specified, can in fact be the desired 
shape of the downstream boundary, in which case the so-called downstream boundary 
becomes the upstream boundary, and its shape is not specified, but depends on the spe- 
cified velocity distributions on the boimdary of the entire flow field. This feature of re- 
versible flow means that, for example, either the upstream or the downstream shape of 
flow area, but not both, can be specified to have constant radii at hub and/or shroud, so 
as to match the axisymmetric surfaces of turbomachine annuli. 

In figure 12, if r^ is zero, a singularity exists at the origin, where 17 takes 
on its full range of values. However, in <^,^,17- space, where the design problem is 
solved, the singularity disappears with the 17 values equally spaced along the boundary 
line for ip equal to zero at constant cp . (Along this boundary line* the velocity is then 
constant . ) 

This polar arrangement of fluid lines on the upstream boundary, with its associated 
singularity, becomes essential when the duct includes a centerbody that starts down- 
stream from the duct inlet. For this case, the stagnation streamline at the upstream 
boundary is located at the origin (r^ = 0.0, in fig. 12), where ip is 0.0 and 17 has 

its full range of values. Downstream from the stagnation point (q.v.) on the center- 
body, the value of ip is then ever3nvhere 0.0, and a grid of 17 and (p lines covers the 
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centerbody surface. This situation is shown in figure 13 for the simple case of incom- 
pressible, potential flow around a sphere. 

Of course, it is not necessary that the upstream boundary contour conform to the 
fluid lines of constant ip and tj in figure 12 . In fact, an arbitrary shape, such as that 

in figure 11, could equally well be treated by this polar grid arrangement, in which 
^ 2 

case, the integrand r^ of equation (30 i) would vary with 0*^ along the line integral. 

Specified velocity distribution on boundary . - The specified velocity distributions 
on the upstream and downstream boundaries are normally uniform, unless either or 
both of these surfaces are planes of symmetry (q.v.), in which case, the velocity is not 
specified, and instead, the boundary condition becomes 9 In q/9c^ =0.0. For those 
cases in which the velocity is uniform, the nondimensional upstream velocity q^ is 
1.0, and the downstream velocity q^^ can have any value, provided it is not super- 
sonic . 

Along the duct wall, the velocity distribution is prescribed as discussed previously 
in the section General Considerations. In addition, for the polar grid arrangement of 
fluid lines on the upstream boundary, if A0jj (fig. 12) is 360®, then, on the two sur- 
faces of constant 17 equal to 0.0 and 1.0, the boundary condition is that the velocity 
(which is not input) is the same at corresponding <77 , z^-points . 

For planar symmetry (q.v.) parallel to the direction of flow, it is necessary to de- 
sign only the part of the duct that lies between planes of symmetry, or in the case of a 
single plane, that lies on one side. A plane of symmetry then becomes a boundary of 
the flow field, and its boundary condition is that the normal gradient of In q is zero. 
For a rectangular grid of fluid lines on the upstream boundary, there are two conven- 
ient planes of symmetry, one of constant ip, the other of constant tj. For the polar 
grid arrangement, the included angle A0* (fig. 12) between two planes of symmetry can 
have any value that divides evenly into 360®. hi these cases of multiplane symmetry, 
the line of intersection between planes of symmetry is the x-axis in x,y, z-space, and, 
only if the flow field does not include this axis, can the duct turn. 

For branched ducts, the flow field in cp,ip,r}- space can be cut by any specified line, 
of constant shape on ?/), 77 -planes, moving directly upstream in the negative (^-direction 
from the downstream boundary. Because the Jp and 17 values on the two surfaces of 
the cut are constant in the direction of <pt these surfaces are composed of streamlines. 
This cut is extended upstream to the value, or values, of (p along which the flow 
branches . The cut line is defined by 






(31a) 

(31b) 
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Vq ~ VQipQ) 



where is path length along the cut line on ^,T]-planes in 9?, 2/^,i7-space. Thus the 
velocity distributions on each of the two surfaces of the cut are prescribed separately 
by 


In =f(9?,p^) ( 31 c) 

hi the region where the flow branches, a stagnation condition (q.v.) can exist, in which 
case, the prescribed velocity is zero. In the absence of stagnation, the prescribed 
velocities remain finite, and the boundary surface at the branch junction is cusped in 
X, y, z-space, Any number of branches for the same duct can be treated by the method 
just described. For any branched duct, the specified velocity distribution should be 
such that, in X, y, z-space, the branch turns away from, and not into, the duct or other 
branches. This condition is achieved by specifying higher velocities on the branch wall 
opposite the duct or other branches. 

For ducts with centerbodies, the polar arrangement of fluid lines on the upstream 
boundary in x,y, z-space is required, as discussed previously. In this case, the origin, 
with ip equal to 0.0 and with 17 having its full range of values (0.0 to 1.0), becomes a 
boundary surface in ^,T]-space. This surface is the stagnation (or cusp) streamline 
upstream from the centerbody and becomes the surface of the centerbody itself down- 
stream from the stagnation point or cusped region. Upstream from the centerbody, the 
boundary condition is constant In q (value not specified, but related by the governing 
differential equation to the values of In q at all adjacent grid points) for all values of rj 
having the same value of (p . Downstream from the stagnation point, or from the cusp, 
the boundary condition is the specified distribution of In q as a function of cp and rj 
on the surface of the centerbody. 

Turning angle . - The turning angle of a three-dimensional duct is given by the 
change in flow direction between the upstream and downstream boundaries . These flow 
directions are given by the direction cosines: cos Q! ^lU* the 

upstream boundary (9? = 0.0) and cos cos and cos at the downstream 

boundary ((p = (p-^). Although, for the design of two-dimensional ducts, analjrtical 
methods have been developed (ref. 1) for predicting the turning angle as a function of 
the prescribed velocity distribution, no attempt to do this for three-dimensional ducts 
is made in this report. 


External Flow Fields 

External flow fields are flow fields that surround bodies in x, y-z-space. The bod- 
ies can be solid or ducted. The three-dimensionality of these flow fields implies that 
the bodies are unsym metrical, although the bodies can have planar symmetry. 
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Upstream boundary . - For external flow fields, the upstream boundary in x,y, z- 
space can be, for convenience, similar to that for those internal flow fields (q.v.) 
where a polar arrangement of fluid lines is used. However, for external flow fields, in 
the absence of planar S5rmmetry, the shape of the upstream boundary in x,y, z- space 
can always be circular with (fig. 12) and r^j 

spectively. 

From equation (30c), with 

yu='^U,max 

^U = ’^U,max 

the nondimensional area of the upstream boundary becomes 

^U = 10=^/4,max0»U <32b) 


L 


equal to 360® and 0.0, re- 


(32 a) 




from which 


U, max 


Vt 


(32c) 


If, by definition, 


»'u 

7J = 

2 tt 


(32d) 


so that the full range of t] values goes from 0.0 to 1.0 as 0^ goes from 0.0 to 27 t, 
then, from equations (30 f) and the continuity equation (9d), 




or, because Py = 1.0, qy = 1.0, and Is 90°, 


_ 1 


d^ = 27rr^ dr^ 


from which. 
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= 27t 


r^v 

Jo 


or 


ij) - irr 


U 


From equation (32e), with r^j given by equation (32c), 




max 


= 1.0 


(32e) 


(32f) 


so that the full range of 7 p values also goes from 0.0 to 1.0, as r^ goes from 0.0 to 

lA/i: 

For the given arrangement of fluid lines, the shape of the upstream boundary with a 
grid of j/j and 17 lines is that shown in figure 14, where the scale of y^ and 2^ is 
such that the upstream area A-^ is 1,0. In this figure, the stagnation (or cusp) stream- 
line is located at the origin, where ijj is 0.0 and 77 has its full range of values from 
0.0 to 1.0. This singularity in x,y, z- space disappears in <73, 77- space, where the de- 
sign problem is solved, and where the 77 values are finitely spaced along the boundary 
for ^ equal to zero. Boundary conditions along this stagnation streamline can be 
treated in a similar fashion to that for the stagnation streamline upstream from the 
centerbody in an internal flow field (q.v.), (Note that the fluid lines of constant 77 in 
figure 14 need not be equally spaced.) Downstream from the stagnation point or cusp 
region on the body, the value of ijj is everywhere 0.0 on the surface of the body, which 
surface is then covered by a network of 77 and cp lines . In this circular grid arrange- 
ment, unlike that for internal flow (fig. 12), the fluid lines of constant if) are spaced 
closer together near the origin, as shown in figure 14. This closer spacing is used be- 
cause the fluid in these stream laminae experiences greater deformations, and, there- 
fore, requires closer attention, in passing around the body. (It is for this reason also 
that the fluid lines of constant 77 might be more closely spaced relative to one side of 
the body.) For external flow fields with this circular grid, it is noted that, for constant 
values of ip, the spacing of interior grid points from the boundary surface is conven- 
iently constant. Although the nondimensional radius r^., of the grid is constant at 

VI/tt, the dimensional value must be sufficiently large, relative to the size of the body, 
that the nondimensional velocity q along the surface of the cylinder in x, y,z-space is 
essentially constant and equal to q^, that is, 1.0. 

Grid configuration . - The boundary of the grid configuration in x, y, z-space is a 
right circular cylinder, as shown in figure 15, with a polar grid at the upstream bound- 
ary, as shown in figure 14. The configuration of the downstream boundary varies with 
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such conditions as closure (q.v.) and planar S 3 anmetry normal to the direction of undis- 
turbed flow. The body itself is located in a small region close to, and including, the axis 

of the cylinder. The grid surfaces of constant (p, ip, and 77 are warped (not shown) in 
the disturbed region of flow around the body. In this disturbed region, (p surfaces 
should be more closely spaced (as shown) for greater detail, as is the spacing of ip 
surfaces in figure 14. 

The boundary of the grid configuration in (p, ?/), 77 - space, which grid corresponds 
point for point with that in x, y, z-space, is a rectangular parallelepiped, as shown in 
figure 15, In this grid, unlike that in x,y, z-space, the grid points themselves occupy 
the corners of small rectangular parallelepipeds. This arrangement results in conven- 
ient finite-difference stars (q.v.) for the solution of the governing differential equation 
in ^,^, 77 - space. 

Specified velocity distribution on body . - The specified velocity distributions on the 
upstream and downstream boundaries are uniform, unless the downstream boimdary is 
a plane of symmetry (q.v.),, in which case, the velocity is not specified on the plane of 
symmetry, and, instead, the boundary condition becomes 9 In q/dcp =0.0. For those 
cases where the velocity is uniform, the nondimensional velocity, both upstream (q^) 
and downstream (Qj^» is 1.0. The dimensional value of this velocity must not be super- 
sonic . 

Provided the nondimensional radius of the essentially cylindrical boundary (fig. 15) 
in x,y, z-space, which radius is approximately equal to r.. (equals >/r 7 ^), is suf- 

w 9 xhcL^ 

ficiently large relative to the size of the body, the velocity on this surface of constant ip 
equal to 1.0 (fig. 14) is uniform and equal to 1.0. If (fig. 12) is 360^, that is, 
there is no planar symmetry parallel to the direction of flow, which is discussed later, 
then the velocity (not specified) is the same at corresponding (^,jp=points on the two sur- 
faces of constant 77 equal to 0.0 and 1 . 0 . 

The remaining boundary is the 77 , (p-surface of constant ip equal to zero, which sur- 
face is the stagnation (or cusp) streamline upstream from the body, and becomes the 
body surface itself downstream from the stagnation point or cusped region. Upstream 
from the body, the boundary condition is that In q (value not specified, but related by the 
governing differential equation to the values of In q at all adjacent grid points) is con- 
stant for all values of 77 having the same value of (p , Downstream from the stagnation 
point or cusp, the boundary condition is the specified distribution of In q as a function 
of q) and 17 on the body surface. And downstream from the body, the boundary condi- 
tion is that In q (value not specified, but see previous discussion for upstream stagna- 
tion streamline) is again constant for all values of 77 having the same value of (p . The 
shape of this boundary downstream from the body in x, y, z-space depends on the speci- 
fied velocity distribution on the body surface. The downstream boundary should degen- 


54 


erate into the downstream stagnation streamline, but whether or not it does involves 
the problem of closure, which is considered briefly later in this report. 

For planar symmetry (q.v.) parallel to the direction of undisturbed flow, that is, 
parallel to q^, only a sector of the cylinder bounded by two planes of constant 0’ in 
X, y, z-space need be designed. These planes of symmetry then become two boundaries 
of the flow field, and their boundary condition is that the gradient of In q with respect 
to 77, that is, normal to the planes, is zero. The included angle A0’ (fig. 12) between 
the planes of symmetry can have any value that divides equally into 360°; in particular, 
for A 6’ equal to 180°, the flow field is divided equally, with each half being a mirror 
image of the other . 

A ducted body results if the flow field in 17- space (fig. 15) is cut along any 
plane surface of constant ipQ greater than zero. This cut starts at the downstream 
boundary and extends upstream in the negative (p-direction to the value, or values, of (p 
corresponding to the upstream stagnation or cusped region of the body. Separate veloc- 
ity distributions on each of the two surfaces of the cut are prescribed by 

In q^ — cp (33) 

As a result of this cut and of the separate velocity distributions, the surface of the duc- 
ted body in x,y, z-space has the constant value of and is covered by a network of rj 
and cp lines. As is the case for unducted bodies, closure (q.v.) depends on the veloc- 
ity distribution prescribed on the body surface. Note that a centerbody exists within 
the duct if a velocity distribution is specified on the centerbody surface (?^ = 0.0), as 
described in the section Specified velocity distribution on boundary, under Internal Flow 
Fields. Otherwise the boundary condition for ip equal to zero is that In q (value not 
specified, but see previous discussions of boundary conditions for stagnation stream- 
lines) is constant for all values of 77 having the same value of cp . 

Closure . - The ip surfaces equal to 0.0 or ipQ define the body shape in x, y, z- 
space. These surfaces should, therefore, close at the downstream end of the body in 
x,y, z-space. Depending on the specified velocity distribution on the body itself, these 
2 p surfaces may close, or they may remain as finite-size (positive or negative) tails of 
the body, extending downstream to infinity. This closure condition for the body is not 
treated in this report. However, for flow fields that are symmetrical about a plane 
normal to the direction of undisturbed flow, which plane is halfway along the length of 
the body, closure is ensured by the symmetry. Also, if the tail is of modest cross sec- 
tion relative to the body size, it can be chopped off with no significant effect on the pre- 
scribed velocity distribution. 
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Numerical Procedure 


The following list is a brief outline of the numerical procedure for the design of 
three-dimensional flow fields. For details on special situations, refer to previous sec- 
tions in this report. 

(1) For internal flow fields, specify the shape of the upstream boundary, that is, 
the upstream flow area of the duct, which boundary is the y, z-plane for x equal to 
zero in x,y, z-space. Compute the area of the upstream boundary, and express the y, z- 
coordinates of this perimeter in nondimensional units, as described in this report. For 
external flow problems, step (1) is not required. 

(2) Prescribe the desired distributions of nondimensional velocity q as functions 
of nondimensional path length s along streamlines (constant ip and tj) on the undefined 
boundary of the flow field in x,y, z-space. These distributions with respect to s are 
easily converted into distributions with respect to velocity potential cp by the equation 


(p = 



q(s)ds 


Thus Q, which is In q, is a known function of cp along streamlines on the lateral 
boundary. The maximum value of specified velocity, when expressed dimensionally, 
should not exceed the local speed of sound (except perhaps slightly, and then only for a 
restricted part of the boundary) . 

(3) Assume initial values for Q at all interior grid points of the finite-difference 

grid. 

(4) Compute the coefficients C from equations (28f), where, for the initial finite- 
difference solution of the governing equation, dn/dip^ dm/dr], and 6 can be assigned 
values of 1.0, 1.0, and 90^, respectively, for a rectangular arrangement of the fluid 
lines on the upstream boundary. (Also, for the initial finite -difference solution, with a 
rectangular arrangement of the upstream fluid lines, the total curvatures and 

are assumed to be zero.) For a polar arrangement of the upstream fluid lines, special 
care is required as the radius r approaches zero. 

(5) At every internal grid point, solve the finite-difference form (eq. (28 e)) of the 
governing differential equation (17h) by using standard finite-difference procedures . 
Equation (28 e) is solved when the highest absolute value of residual error (R in the flow 
field is less than a small specified amount. 

(6) At each interior grid point of the finite-difference grid, compute the direction 
cosines from equations (21c), (21d), (2le), (22c), (22d), (22e), (23c), (23d), and (23e). 

(7) Compute the total curvatures and from equations (17 f) and (17 g), re- 
spectively, or from equations (C5) and (C6), respectively. 

(8) Compute the angles 6 from equation (8c). 
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(9) Determine the densities p from the velocities q. 

(10) Determine the physical shape of the flow field in x,y, z -space by computing the 
x,y, z -coordinates for every point in the finite -difference grid by using equations (26a), 
(26d), and/or (26g) for x, equations (26b), (26e), and/or (26h) for y, and equations 
(26c), (26f), and/or (26i) for z. 

(11) Compute the parameters A and B from equations (24) and (25), respectively. 
Although in theory the values of A and B are determined by equations (24) and (25), 
respectively, in practice, the values may be better determined, during the iteration 
process, by equation (lOd) and the ratio of equations (24) and (25) or by equation (lOd) 
and the difference of equations (16c) and (16d), which difference gives 



(12) Repeat steps (4) to (11) until the value of velocity at each interior grid point has 
converged within a specified limit or until the maximum absolute value of the residual 
(ft for a new set of coefficients C from step (4) is less than a small specified amount. 


CONCLUDING REMARKS 

The general design method developed in this report is for steady, three-dimen- 
sional, potential, incompressible or subsonic -compressible flow, hi this design meth- 
od, the flow field, including the shape of its boundary, is determined for arbitrarily 
specified, continuous distributions of velocity as a function of arc length along the 
boundary streamlines . The method applies to the design of both internal and external 
flow fields, including, in both cases, fields with planar symmetry. These designs re- 
sult from the finite-difference solution of a governing differential equation for the dis- 
tribution of velocity in transformed space, the coordinates of which are the velocity po- 
tential and two steam functions . The anal 3 d;ic problems associated with stagnation 
points, closure of bodies in external flow fields, and prediction of turning angles in 
three-dimensional ducts are discussed, but not treated in detail. 

This three-dimensional design method applies to simple ducts, branched ducts, an- 
nular ducts, ducts with centerbodies, including the nose and/or tail region, and ducted 
or nonducted bodies in infinite, uniform flow fields. Solutions determine a curvilinear 
coordinate system in physical space for boundaries with prescribed and, therefore, 
basically good velocity distributions, which coordinate system, if desired, can then be 
used for subsequent, three-dimensional, viscous flow analyses. 
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APPENDIX A 


ORTHOGONALITY OF STREAM SURFACES AND POTENTIAL SURFACES 

Two families of surfaces are orthogonal, if the dot product of the gradients of their 
generating functions is zero. Thus, for potential surfaces of constant <p and stream 
surfaces of constant ip. 


\ 8x 9y 9z / \ 8x 9y 9z 


9x 9x 9y 9y 9z 9z 


or, from equations (5d) and (6c), 


vcp ■ ^ M + /^ ^ - M mW + /m. 

y9y 9z 9z 9yy 9x l9z 9x dx dzjdy y9x 9y 9y 9x/ 9z 


= 0 

so that stream surfaces of constant ip and potential surfaces of constant cp are or- 
thogonal. 

Likewise, for stream surfaces of constant 7) and potential surfaces of constant (p, 

V„ . Vcp + 

y 9x 9y 9z y y 9x 9y 9zy 

dx dx 9y 9y 9z 9z 
or, again from equations (5d) and (6c), 

p vr) ■ vcp ^ - M + ^ ^ 

9x y9y 9z 9z 9yy 9y \9z 9x 9x dzj 9z \9x 9y 9y 9x 

= 0 

so that stream surfaces of constant 77 and potential surfaces of constant (p are or- 
thogonal . 
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For the two families of stream surfaces 


dx dx dy dy dz 9z 


which, in general, is not zero, so that the stream surfaces of constant ip and the 
stream surfaces of constant 77 are not orthogonal. In fact, from equations (lOd), (lOf), 
and (lOh), 


Vip • ^ 

tan 0 


(Al) 


so that stream surfaces of constant ip are orthogonal to stream surfaces of constant rj 
only when 0 is 90° (a truism). 

However, because both families of stream surfaces are normal to the potential sur- 
faces of constant (p, the angle G (fig. 3) with which the stream surfaces intersect is 
measured on the potential surface and is given by 


cos 6((p,ip,r}) =— 

ysG 

where 



9x , 9y 9y ^ dz 9z 
dip drj dip dj] dip dr] 



The angle 6(cp,ip,ri) is also given by equation ( 8 c). 


(A2) 


(A3) 
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APPENDIX B 


GEODESIC CURVATURES y OF (p, AND rj CURVILINEAR COORDINATE 
LINES ON jp AND r] STREAM SURFACES IN x,y, z-SPACE 

The local geodesic curvature of a coordinate line on a stream surface measures 
the local deviation of that line from the shortest distance on the surface between two 
local points an elemental distance apart on the line. (On a plane surface, the local geo- 
desic curvature of a line is the reciprocal of its local radius of curvature,) From page 
295 of reference 6, on a surface with u, v curvilinear coordinates, a field of curves 
having the imit tangent vector e(u,v) has at every point the geodesic curvature 


y = n ■ Vg X e 


i.1 

H 




8v 




(Bl) 


where n is the unit vector normal to the surface, r is the vector position on the sur- 
face, Vg is the vector differential operator on the surface, and 


H =,§£x-^ 
9u 9v 


n 


(B2) 


First, consider the t] 


stream surface shown in sketch (b). 



For this surface. 


- - - 
n = e^ X eg 

u = <p ) 


(B3) 


V = Ip 


J 
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From equations (lOi), 


8u d(p d<p ^ q ^ 


_ 1 — 

av dip dip ^ qB ^ 




y 


(B4) 


Vq = v</} -^ + 'V^ip — 
^ d(p ^ dip 


yds -y d(p 1 dn dip 


— a , -Q— a 

= qe_ — + qBe^ — 
^ dcp ^ dip 


(B5) 


and, from equations (B2), (B3), and (B4), 


H = 


A e \ X e \ • ^ X 62) = — 

V / q"B 


(B6) 


Thus, for the (^-coordinate curve on the 77 stream surface, with e equal to e^, equa- 
tions (Bl), (B3), (B4), (B5), and (B6) give 


de __ ae . 2 

7l = (e^ X e^) • (qe^ ^ ^ ® 




which reduces to 


Il = e. 


ae 


_ g a In q 
d(p dip 


(B7) 


so that 


dcp dip 


ase) 
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Equation (13e) is the irrotationality condition on the tj stream surface. 

Likewise, for the ^-coordinate curve on the ?]-surface, with e equal to 62 * 
equations (Bl), (B3), (B4), (B5), and (B6) give 

.Zl = Be ^^2 _ e In q , 8 In B 
q ^ dip dtp dtp 


However, 


so that 


and, thus. 


or 


e 


1 


= 0 



l2^g- ^^1 . 3 In q 6 In B 

q ^ dip dtp dtp 



^ a In q 3 In B 

dip dtp dtp 


Equation (16c) is the continuity condition for flow in j] stream laminae. 
Next, consider the ip stream surface shown in sketch (c) . 


(B8) 


(16c) 
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For this surface 


From equations (lOi), 


n = e^ X ^ 


U = T) 

V = cp J 


> 


(B9) 


0r ^ 

= ^ = 

dm — 

■ ■ ft 

= J-e 

8u 

'dr} 

dr, 3 

qA = 

dr 


ds 

= ' ft 

= ie 

dv 

dtp 

dv 1 

q ^ 




(BIO) 




\ dm j dr} yds dcp 


— 9 ^ 9 

= Qi^^e^ — + qe^^ — 


dr} 


dtp 


(Bll) 


and, from equations (B2), (B9), and (BlO), 


H = 


q^A 


(B12) 


Thus, for the 77 -coordinate curve on the Jp stream surface, with e equal to e^, equa- 
tions (Bl), (B9), (BIO), (Bll), and (B12) give 


^3 _ _ 9 In q 8 In A 

q ^ dr} dtp dtp 


8e, 


(Bl^ 


or 
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Ae . - 8 In q 3 In A 

^ dri Bcp d(p 

Equation (L6d) is the continuity condition for flow in ^ stream laminae . 

Likewise, for the ^-coordinate curve on the ^-surface, with e equal to e 
tions (Bl), (B9), (BlO), (Bll), and (B12) give 

q ^ dcp drj 


or 

9t) dcp 

Equation (13f) is the irrotationality condition on the ip stream surface. 


aed) 
, equa- 


(B14) 


asf) 
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APPENDIX C 


TOTAL CURVATURES K OF ip AND 7] STREAM SURFACES IN x,y, z-SPACE 

From page 306 of reference 6, the local total curvature of the tj stream sur- 
face is given by 

. r\ 

^ / dip\d<p ) 

where 

R=Yie^ +Y2«2 

and H, 9r/9j^, 9r/9^, and are given by equations (B6), (B4), (B4), (B7), and 
(B8), respectively. Thus 


(Cl) 


(C2) 



which reduces to 


K /9e„ 

q2 




a7f) 


Equation (L7f) gives the value of K^/q in the governing differential equation (17h). 
Likewise the total curvature of the ip stream surface is given by 


where 





(C3) 


(C4) 


and H, dr/dcp, 9r/9rj, y^* 7^ are given by equations (B12), (BIO), (BIO), (B13), 

and (B14), respectively. Thus 
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which reduces to 


S ^ 


dn 




_ 9 _ 

dcp 




K, 

-- i= A 


2 

q 



(I7g) 


2 

Equation (17 g) gives the value of K^/q in the governing differential equation (17 h) . 

From equations (B7), (B8), (B13), and (B14), the geodesic curvatures y, and, 
therefore, the total curvatures K, can be expressed in terms of velocity rather than 
unit vectors. Thus, for the t] stream surface, from equations (Cl), (C2), (B4), (B6), 
(B7), and (B8), 




which expands to 


^77 8 In B/a In q . 

q2 dcp \ dcp 


d In B\ 
dcp ) 


d^ Inq 
dcp^ 


d^ In B ^2 


9 In B 8 In q ^ 

dip dip 


d^ In q\ 

dip^ J 


(C5) 


Likewise, for the ^ stream surface, from equations (C3), (C4), (BIO), (B12), (Bl3), 
and (B14), 



which expands to 


^ij) 8 In A /8 In g I 8 In A \ 9^ In q 9^ In A ^2 1 3 In A 3 In g , 8^ In q\ 
q2 dcp \ dcp dcp / 9^2 g^2 ^ 9^ dr; ^2 J 

If equations (C5) and (C6) are introduced into equation (17h), the governing differential 
equation becomes 
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8^ In q , 9^ In p , 8^ In sin 6 8^ In B 8^ In A 

dcp^ d(p^ dcp^ dcp^ dcp^ 

from which, assuming the constants of integration are zero, 

p sin 9 = qAB 

which is the continuity equation (lOd) . 



APPENDIX D 


ALTERNATIVE DERIVATION OF GOVERNING 


PARTIAL-DIFFERENTIAL EQUATION 


In addition to its derivation in the main text of this report, the governing, second- 
order, partial- differential equation for the distribution of In q in cp , if), r]- space can be 
derived from separate considerations of the two-dimensional continuity and irrotation- 
ality conditions on the t) and ip stream surfaces, the results of which are then inter- 
related by the simple continuity equation (lOd) for three-dimensional flow. 

From the two-dimensional, continuity equations (16c) and (16d), or, from equa- 
tions (B8) and (B13), which are the same. 


3 In B /a In o , 3 In B\ , p ^^2 ^^1 , ^ 3^ In g 3^ In B 

3(p \ 3(p 3<p / 8(p dip ^ dip dtp g^2 g^2 

_ 2 

6 In A /e In q , 3 In AV ^ _ 0^ In q 8^ In A 

d(p \ d(p dip ) dip dq ^ dip dq 


pi) 


(D2) 


Equations (Dl) and p2) are the two-dimensional, continuity equations in the q and ip 
stream laminae, respectively. 

From the two-dimensional, irrotationality equations (13e) and (13f), or, from equa- 
tions (B7) and (B14), which are the same. 


9e„ de 


Be. 


;\ 2 — 
o e 


dip dip 


dip dip 


1 _ j^2[ 8 In B 8 In q , 8^ In q 


dip dip 


P3) 


dip‘ 


dq dip 


+ Ae, 


2 — 

^ ®1_ ^ ^2 /a In A 81n q ^ 8^ In q ] 


dip dq 


dq dq 


m 


dq^ 


Equations (D3) and (D4) are the two-dimensional, irrotationality equations on the q 
and ip stream surfaces, respectively. 

Combining equations pi) and pS) to eliminate d^j^/dip dip and introducing the 
total curvature from equation (17f) give 
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(D5) 


a In B /a In q ^ a In b\ ^ ^ ^ jj2 [ 8 In B 8 In q ^ 8 ^ In q \ ^ 8 ^ In q ^ 8 ^ In B 

8 (/? \ 8 <p dcp ) 2 \ dip dip ^2 j g^2 g^2 


Equation (D5) is the equation for two-dimensional flow in 7 ] stream laminae . 

O 

Combining equations (D2) and (D4) to eliminate 8 e^/^cp 817 and introducing the 
total curvature from equation (17 g) give 


8 In A /a In q , 8 In A\ _ 5^ ^ ^2 p In A 8 In q , 8^ In q \ _ 8^ In q 8^ In A 

8<^ V 8(p dcp ) 2 y &n ^ 9^2 y 9^2 9^2 


(D 6 ) 


Equation (D 6 ) is the equation for two-dimensional flow in ip stream laminae. 
From the simple continuity equation (lOd) for three-dimensional flow. 


8 ^ In q 8 ^ In A 8 ^ In B _ 8 ^ In p 8 ^ In sin 6 

dcp^ dcp^ dcp^ dcp^ dcp^ 

from which, after equations (D5) and (D 6 ) are added, 


8 ^ In q ^ 3 ^ In p ^ 8 ^ In sin 0 
dcp ^ dcp ^ dcp ^ 


K. K 


Ijk - ^ + ^2( 8 In B 8 In q , 8 " In q 


dip 


dip 


dlp‘ 


+ A 


2 


8 In A 8 In q ^ 
dr] dr] 


8 ^ In q\ 

877^ / 


8 In B /s In q ^ 8 In B 
dcp \ dcp dcp / 


8 In A /8 In q ^ 8 In a \ _ ^ 
dcp \ dcp dcp / 


(17h) 

Equation (17h) is the governing differential equation for the distribution of In q in 
( 7 ?, ^,77~space as developed previously in the main text of this report. 


69 



REFERENCES 


1. Stanitz, John D.: Design of Two-Dimensional Channels with Prescribed Velocity 

Distributions Along the Walls . NACA Rep. 1115, 1953. (Supersedes NACA TN 
2593 and TN 2595.) 

2. Yang, Tah-teh; Hudson, William G.j and Nelson, Carl D. ; Design and Experi- 

mental Performance of Short Curved Wall Diffusers with Axial Symmetry Utilizing 
Slot Suction. NASA CR-2209, 1973. 

3. Yih, Chia-Shim: Stream Functions in Three-Dimensional Flows. La Houille 

Blanche, vol. 12, no. 3, July- Aug. 1957, pp. 439-444. 

4. Milne- Thomson, L. M.: Theoretical Hydrodynamics. The Macmillan Co . , 1950. 

5. Phillips, H. B.: Vector Analysis. John Wiley & Sons, Inc., 1933. 

6. Brand, Louis; Vector and Tensor Analysis . John Wiley & Sons, Inc., 1947. 


70 



Potential surface with 



71 



ijl + (tjl 



72 


,-igure 3. Structure of flow field in x, y, z-space, showing velocity 
vector q;(^, Tj potential surface of constant (pitp.ip and rj.tp 
stream surfaces of constant r\ and q/, respectively; vector gra- 
dients of (p, iji, and V; angle 6; and differential path lengths 
ds, dn, and dm. 















Figure 5. - Sketch of vector gradients 7tp and Trj of stream functions on potential surface 
of constant (p, showing relations among these gradients, the angle 0, and various unit 
vectors. 
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7} + d7J 



Figure 8. - Sections of three-dimensional stream laminae resulting from inter- 
sections of adjacent stream surfaces of constant $ and t? with potential sur- 
face of constant <p in this figure; velocity q is normal to these sections. 



Figure 9. - Shape of flow field in transformed ^.(^.rj-space, showing 
spacings ai, 32, and between adjacent points in finite-difference 
grid formed by surfaces of constant (p, tp, and i?. 
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Figure 10. - Finite-difference star in transformed ip,$,Ti-space 
with central grid point 0 and adjacent grid points 1 to 6; vari- 
able spacing between adjacent grid points. 



Figure 11. - Arbitrarily specified shape of upstream boundary 
for internal flow field in x, y, z-space. For nondimensional 
variables, as defined in report, ip and v equal yy and 
Zy, respectively, when 0 is 90° on upstream boundary. 
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Figure 14. - Polar arrangement of fluid lines on upstream boundary for external 
flow fields in x,y,z-space. 
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ment of fluid lines on upstream boundary in physical flow field. 
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